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Capacitor  for  direct  oscillator  frequency  control) . 
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Item  #  19  continued 

quencles  of  the  SC-cut  resonators  Indicate  that,  over  an  operational  temperature  range 
of  -55°C  to  -f-lOO^C,  a  clock  accuracy  of  several  mlllisecond/day  should  be  attainable; 
higher  accuracy  probably  feasible  with  further  development  and  greater  sophistication. 

2.  Severe  B-mode  activity  dip  at  some  temperature  within  the  test  range  of  -55°C  to 
+100°C  noted  In  most,  but  not  all  SC-cut  crystals  used  In  the  test  program. 

3.  Simple  compensation  for  thermal  shock  effects  shown  to  be  beneficial  and  essential 

to  meeting  above  accuracy  goal;  practicality  of  compensation  for  residual  hysteresis  effects 
will  require  further  study. 

4.  Partial  inadequacy  of  polynomial  approximation  for  compensation  purposes;  If  ultimate 
accuracy  Is  to  be  achieved,  operational  ECXO  will  probably  require  auxilary  look-up  table 
In  combination  with  a  low  order  polynomial  approxioiatlon  algorithm. 

5.  Compensation  capability  does  not  appear  to  deteriorate  significantly  In  low  temperature 
region  (in  spite  of  the  severe  C-mode  frequency/temperature  slope  at  the  low  temperature 
extreme) . 

6.  The  phase-modulated  resonance  reflectometer  used  in  this  investigation  offers  a  new, 
valuable  Instrumentation  technique  In  present  and  future  crystal  testing  applications. 
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SECTION  1 


PROGRAM  DESCRIPTION 

1.1  PROGRAM  BACKGROUND  AND  OBJECTIVES 

Advanced  military  communication ,  navigation,  and 
identification  systems  are  placing  increasing  demands  for  higher 
accuracy,  lower  power  dissipation,  smaller  size,  and  greater 
reliability  in  component  timing  and  frequency  control  systems. 

The  unique  problem  of  attainable  clock  accuracy,  for  example,  can 
play  a  deciding  role  in  the  specification  and  selection  of  any 
newly  proposed  identification  system  or  communication  system 
using  spread  spectrum  techniques.  Any  proposed  tactical  system 
based  on  a  clock  accuracy  that  can  only  be  obtained  with  an 
atomic  frequency  standard  is  generally  unacceptable.  Similarly, 
the  power  requirements  of  the  conventional  oven-controlled 
crystal  oscillator  may  be  excessive  in  battery-operated 
applications . 

More  than  two  decades  of  development  effort  have  been 

directed  toward  the  development  of  temperature-compensated 

1  2 

crystal  oscillators  (TCXO)  '  .  These  ovenless  oscillators 
attempt  to  maintain  frequency  accuracy  over  a  moderate 


(1)  O.J.  Baltzer,  "Temperature  Compensation  of  Transitor ized 
Crystal  Oscillators",  12th  Southwestern  IRE  Conference,  Houston, 
Texas;  April  22,  1960. 

(2)  V.  Rosati,  S.  Schodowski,  and  R.  Filler,  "Temperature 
Compensated  Crystal  Oscillator  Survey  and  Test  Results",  Proc . 
37th  Annual  Preouencv  Control  Symposium;  1983,  pp.  501-505 


temperature  range  by  using  a  temperature-sensing  element  and  a 
voltage-variable  capacitor  (varactor)  to  shift  or  ”pull"  the 
oscillator  frequency  in  a  direction  to  counteract  the  natural 
frequency  drift  of  the  crystal.  Currently  available  TCXO 
oscillators  can  employ  analog,  digital,  or  even  microprocessor 
elements  for  the  compensation  circuitry.  The  basic  accuracy  and 
stability  of  all  conventional  TCXO  devices  is  limited,  however, 
by  the  varactor  and  the  basic  TCXO  concept  itself. 

This  traditional  approach — forced  internal  tuning  of  the 
oscillator — introduces  operational  and  performance  problems  into 
any  tactical  system  requiring  highest  accuracy  TCXO  timekeeping 
capability.  Fortunately,  there  is  growing  recognition  of  an 
alternative  technique-one  in  which  the  required  tunability  and 
time  adjustment  functions  are  obtained  by  means  external  to  the 
oscillator:  (3) 


fy^curnntly.  She  tyttm  dttigntT  is  sanfyvntsd  'jizh  a  rssd  rc 
oi'ovids  for  frsqumnay  or  tims  (phass)  adjusvabiliiy.  In  zhs  oast, 
mors  oftsn  than  not  it  appsarsd  sasisst  to  ths  systsm  dssignsr  zo 
sps^fy  an  adjustdbls  frsqusnay  standard  or  aleak.  Irystal 
dsviess  can  os  mads  tsarabls  by  adding  a  aapasitor  in  oarallsl  or 
tn  MTiss  to  ths  quarts  orystal  rssonator.'  Tunability  is  than 
aahisv^  by  varying  ths  talus  of  this  aapacitor  msahaniaalZy  or 
slsetriaally  ( varaator) . . , 


In  most  oasst,  hausvsr,  it  is  ovsrlaoksd  t^uxt  ths  addition  of 
tunability^  affs^s  ths  basic  opsratian  of  zns  f'^cusney  stcandarz 
or  alack  in  a  astrimsntal  ^oay •  Irsqusncu  standards  dsrivs  thsir 
high  stability  and  aaeuraay  from  ths  fact  that  ths  ssssnzial ' 
aontrol  slsmsnt,  ^s  quarts  crystal  rssonator  or  idn,s  atomic 
^senator,  has  a  nigh  rssonanas  5,  ad,  that  thsir  ^scnancs 
frsqusnay  is  highly  invariant  ’jith  tins  or  sstsrnal  taramstsr 
anangss.  Ihs  additi’on  of  a  tuning  aasaoitor,  in  -r.s  toss 
crystal  oseillators,  or  rne  provision  for  magnstic  fisld  -yiriabil- 
...jf  in  ths  aass  of  atomic  ""ssonators,  virtvally  aluaus  zso’KLdss 
ths  psrformar.es  by  providiTig  a  dirsat  couoling'  of  ^s  ■^sorator  tc 
•jorying  sstarral  infZusnass . . .  . 


It  ii,  thai^rort,  uatfuL  to  assart  that  Largs  tUTrability  std 
stata—cf^ths—art  par^ jirvBies  in.  fracuaney  stcatdards  std  tlacks  ora 
vnooripaatbLi f  aiui  to  ths  bast  ci  tha  authors '  hnouLadaa  bava  ravar 
baan  aatiafaetorily  aambinad.  In  othar  wo'ris,  va  hats  tha  asia^n 
that  tunatility  oauaas  datarioration  of  slock,  and.  oaeilhxsor 
parformnaa.  Inataad  of  raquirinq  tha  oLoek  rcatufaaturar  to 
supply  a  ttovibl4,  suparpraeiaion  alookf  ayatam  dasipnars  mould 
aonsidtr  implamgnting  timability  by  mans  astamal  to  tha  actual 
alockt  To  this  and,  fracuaney  tunacHity  aoi  ba  aahiavad  zu 
audding  an  astamal  diraet  synthasrizar  or  a  saeond  tunabla 
oaeillator  ’jith  a  aynthaaizing  loop.  If  phasa  or  tuna  adfuatmant 
ia  daatrad,  astamal  phaaa  shifting  by  digital  or  analog  moans 
should  ba  tha  method  of  ahoiea.'* 


This  concept  of  'external  phase  shifting  by  digital  or 
analog  means”  can  obviously  be  extended  to  an  oscillator 
requiring  temperature  compensation.  Should  the  resultant  device 
be  called  a  "TCXO”?  Possibly.  However r  to  avoid  confusion  with 
the  traditional  TCXO  involving  internal  control  by  means  of  a 
varactor >  the  terminology  "externally  compensated  crystal 
oscillator  (ECXO)”  appears  preferable  for  any  class  of  crystal 
oscillator  in  which  frequency  tunability  of  phase/time  adjustment 
is  obtained  by  means  external  to  the  crystal  oscillator. 


(3)  H.  Hellwig;  S.  Stein  and  F.  Wells,  NBS;  A.  Rahan,  RADC 
Hanscom  AFB,  'Relationship  Between  the  Performance  of 
Time/Frequencv  Standards  and  Navigational/Communication  Svstems," 
Anaaal  eiTl,  1978,  pp.  37-53. 


In  late  1979  Tracer  was  awarded  an  Air  Force  contract 
(F41608-78-G-0063  SG01C3)  to  study  the  "Constancy  of  the  Center 
Frequency”  for  the  Onited  States'  concept  of  the  NATO 
Identification  System.  It  soon  became  apparent  that  accurate 
clock  performance  was  a  critical  factor  in  determining  the 
minimum  useable  code  validity  interval  and  related  sychronization 
requirements  in  the  program.  In  search  of  a  suitable  clock 
candidate  capable  of  accurate »  ovenless  operation  over  a  -55°C  to 
+  85°C  range#  Tracor  recognized  the  unique  and  important 
advantages  of  the  ECXO  approach  over  the  traditional  TCXO  design. 

The  earliest  ECXO  design  considered  by  the  Tracor  team  was 
fairly  conventional.  It  was  based  on  an  AT-cut  crystal 
oscillator  with  a  standard  thermistor  used  as  the  crystal's 
temperature  sensing  element.  However#  whereas  in  conventional 
TCXO  units  the  thermistor  is  used  to  control  a  varactor  within 
the  oscillator's  tuning  circuit#  the  Tracor  design  proposed  the 
utilization  of  the  thermistor  to  control  the  clock's  timing 
output  by  means  of  a  digital  divider/phase  shifter  external  to 
the  crystal  oscillator. 

During  this  study#  Tracor  also  became  aware  of  new 
developments  in  SC-cut  crystal  technology#  including  the 
significant#  but  at  that  time#  little  recognized,  exper ir.iental 
work  at  Hewlett-Packard  into  the  use  of  the  temperature-sensitive 


B-node  frequency  of  an  SC-cut  crystal — using  the  ECXO 

concept— for  temperature  compensation  of  the  concurrent  C-mode 
4  5 

frequency.  ' 

Although  this  original  Hewlett-Packard  experimental  study 
was  limited  to  several  laboratory  oscillators  using  the  dual 
B-mode  and  C-mode  of  operation,  with  testing  restricted  to  the 
0°C  to  40°C  range,  the  Tracor  staff  recognized  the  tremendous 
potential  of  this  work.  Accordingly,  in  submitting  a  Final 
Report  to  the  O.S.  Air  Force  on  the  proposed  NIS  tactical 
identification  system,  Tracor  strongly  recommended  that  further 
research  and  development  be  directed  toward  the  ECXO  concept  as 
an  alternative  to  the  conventional  TCXO. 

In  June,  1981,  the  Electronic  Systems  Division,  Air  Force 
Systems  Command,  Hanscom  AFB,  awarded  Tracor  a  contract  (Contract 
P19628-81-C-0087)  to  test  and  evaluate  an  externally  compensated 
oscillator  (ECXO)  utilizing  dual  mode  frequencies  of  an  SC-cut 
resonator  as  the  means  of  temperature  compensation.  The  first 
phase  of  the  contract  program  included  a  study  and  investigation 
of  the  behavior  and  inter-relationship  of  the  dual  B-mcde  and 
C-mode  frequencies  of  various  SC-cut  resonators  under  varying 
temperature  conditions  (with  data  to  be  taken  at  temperatures 

(4)  J.  Kusters,.M.  Fischer,  and  J.  Leach,  "Dual  Mode  Operation  of 
Temperature  and  Stress  Compensation  Crystals",  Proc.  32nd  Annual 
Symposium  on  Frequency  Control,  pp  389-397,  1978 

(5)  U.S.  Patent  4,079,280;  "Quartz  Resonator  Cut  to  Compensate 
for  Static  and  Dynamic  Thermal  Transients",  J.  Kusters,  J.  Leach, 
M.  Fischer;  March  14,  1978 


ranging  from  -55°C  to  +100°C)/  together  v/ith  the  delivery  of  two 
ECXO  test  oscillator  units.  This  research  program  was  performed 
during  the  period  1  June  1981  -  June  1982.  Subsequently,  under 
an  optional  task  statement,  the  contract  was  extended  to  include 
further  investigation  into  hysteresis,  re-traceability, 
stability,  and  aging  characteristics  of  dual  mode  operation  under 
controlled  temperature  cycling  and  intermittent  operation.  This 
portion  of  the  contract  work  covers  the  period  1  January  1983  - 
31  December  1983. 

1.2  ACKNOWLEDGEMENT 

This  work  was  performed  under  the  supervision  of  Dr.  O.J. 
Baltzer,  Principal  Investigator.  Mr.  Charles  S.  Stone, 

Brightline  Corporation,  under  a  consulting  agreement  with  Tracor, 
Inc.,  had  major  technical  responsibility  throughout  the  project, 
including  the  design  of  the  instrumentation  system  used  in  the 
investigation  and  the  design,  subsequent  modification, 
calibration  and  testing  of  the  ECXO  units.  Charles  Baltzer 
developed  the  digital  portions  of  the  ECXO  design,  including  the 
preparation  of  the  compensation  look-up  table  and  syntonizaticn 
algorithms.  Robert  L.  Nelson  was  responsible  for  a  study 
exploring  the  feasibility  of  the  hardware  implementation  of  an 
ECXO  system,  as  an  alternative  to  a  microprocessor 
implementation.  Clement  Ip  undertook  a  large  share  of  tlje 
environmental  testing  and  measurement  of  the  B-mode  and  C-moce 
frequency  characteristics  of  the  various  SC-cut  quartz  resonators 
tested;  he  also  participated  in  the  software  development  of 


various  data  logging  and  data  analyses  programs.  We  also  wish  to 
acknowledge  that  several  of  the  SC-cut  resonators  used  in  this 
project  were  kindly  supplied  by  Dr.  John  Vig,  EPADCOM,  Ft. 
Monmouth,  N.J.  The  authors  are  also  deeply  indepted  to  Mr. 
William  F.  Donnell,  Tracer  Program  Manager,  and  Mr.  Alfred  Kahan, 
RADC/ESE,  for  their  many  helpful  discussions  and  their  assistance 
and  constant  support  throughout  the  contract  program. 
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SECTION  2 


DISCUSSION  OF  THE  ECXO  CONCEPT 
2.1  ADVANTAGES  OP  THE  ECXO  APPROACH 

Virtually  all  of  the  crystal  oscillators  currently  used  in 
communication  applications  attempt  to  achieve  accuracy  in 
time/frequency  by  means  of  internal  frequency  control  ('tuning') 
of  the  oscillator  itself.  In  the  case  of  oven-controlled  crystal 
oscillators  (CCXO) ,  the  quartz  resonator  is  maintained  at  a 
constant  temperature  close  to  its  frequency  turning  point; 
consequently/  only  vernier  frequency  tuning  is  then  required. 
However/  in  the  case  of  ovenless  temperature-compensated  crystal 
oscillators  (TCXO)  utilizing  internal  frequency  control,  the 
problem  is  more  severe.  The  quartz  crystal  is  deliberately 
"pulled"  off  of  its  natural  frequency  resonance,  often  by  an 
extreme  amount,  to  compensate  for  the  temperature-dependent 
characteristics  of  the  quartz  crystal.  The  required  'tunability' 
is  only  obtained  at  the  expense  of  instability  and  loss  of 
accuracy.  As  a  result  of  this  intrinsic  limitation,  even  though 
much  research  and  development  effort  has  been  given  to  the 
conventional  TCXO  approach  for  more  than  two  decades,  currently 
available  TCXO  units  fall  far  short  of  providing  reliable 
frequency  accuracy  over  any  extended  temperature  range  or 
interval  of  time. 

The  external  compensated  crystal  oscillator  (ECXO)  approach, 
in  contrast  to  the  conventional  TCXO,  retains  the  intrinsic 
stability  and  accuracy  inherent  in  quartz  crystal  itself.  This 


is  achieved  by  utilizing  external  adjustment  of  the  output  signal 
rather  than  internal  'pulling'  of  the  quartz  crystal  element. 

Tracer  has  applied  this  principle  (external  compensation  of 
a  crystal  oscillator »  under  microprocessor  control)  for  many 
years  in  various  models  of  the  Tracer  airborne  Omega  navigation 
system.  An  oven  controlled  precision  OCXO  is  used  as  the  clock 
for  the  computer  and  as  the  reference  phase  source  for 
measurement  of  all  received  signal  phases.  The  OCXO  accuracy  is, 
however,  orders  of  magnitude  away  from  that  needed  for  OMEGA/VLF 
navigation.  To  compensate  for  oscillator  inaccuracy,  the 
microcomputer  makes  repeated  measurements  of  the  clock  offset, 
based  on  all  of  the  received-signal  phases.  These  measurements 
are  used  both  to  correct  the  phase  ('time')  of  the  clock  and  to 
determine  its  phase-rate  ('frequency  drift')  error  in  a 
second-order  tracking  loop.  The  corrections  thus  obtained  are 
not  applied  directly  to  the  oscillator,  which  is  allowed  to  run 
freely;  rather,  they  are  employed  entirely  within  the  computer  to 
produce  the  effect,  external  to  the  oscillator,  of  much  greater 
accuracy.  The  present  project  differs  in  part  from  the  CMEGA/VLF 
compensated  OCXO  oscillator  in  that  it  involves  direct 
temperature  compensation  of  an  ovenless  oscillator,  operating 
over  a  wide  temperature  range. 

The  current  feasibility  of  high  accuracy,  low  power, 
ovenless  ECXO  is  largely  made  possible  by  technological  advances 
and  development  in  several  distinct  areas; 

1.  Development  and  availability  of  new  quartz  resonators  (e.c.. 


SC-cut)  with  radically  improved  performance  features; 

2.  Recognition  of  the  B-mode  frequency  (in  SC-cut  crystals)  as  a 
highly  precise,  digital  thermal  sensor; 

3.  Recognition  of  the  external  compensation  approach  as  the 
proper  means  of  achieving  improved  time/frequency  accuracy. 

4.  Significant  advances  in  microprocessor  technology. 

The  properties  and  advantages  of  the  doubly-rotated  SC-cut 
quartz  crystal  are  well  known  and  have  been  widely  described  in 
various  papers  in  recent  years.  Only  a  brief  summary  need  be 
given  here.  The  major  advantages  cited  for  C-mode  operation  with 
SC-cut  crystals,  relative  to  an  AT-cut  crystal,  include: 
improved  frequency-temperature  characteristics  under  both  static 
and  dynamic  conditions;  insensitivity  to  thermal  transients; 
reduced  sensitivity  to  acceleration  (attitude,  shock,  and 
vibration);  and  freedom  from  activity  dips. 

One  potential  problem  with  the  SC-cut  crystal  is  the  E— mode 
of  oscillation,  some  10%  higher  in  frequency  than  the  C-node.  In 
most  oscillator  applications,  the  B-mode  is  unwanted  and 
suppressed — either  by  greater  oscillator  circuit  complexity  or  by- 
use  of  a  crystal  design  which  increases  the  resistance  of  the 
B-mode  of  oscillation  considerably  above  that  of  the  desirec 
C-mode.  In  the  ECXO  application,  however,  the  B-mode  frequenci 
is  purposefully  used  as  a  precise  means  of  measuring  crystal 
temperature  directly. 

Other  disadvantages  cited  for  the  SC-cut~in  non-ECXO 
applications — are  the  manufacturing  difficulties  of  ccntrollirc 


the  temperature  turning  point,  inflection  temperature,  and  the 

finishing  frequency  tolerance  to  the  required  accuracy.  Another 

disadvantage,  in  non-ECXO  applications,  is  the  "stiffness"  of 

SC-cut  crystals,  particularly  of  3rd  or  5th  overtone  types;  it  is 

difficult  to  trim  these  resonators  by  the  conventional  capacitor 

tuning  adjustment  method.  However,  the  compensation  techniques 

used  within  the  ECXO  provides  a  versatile  means  of  frequency 

adjustment,  even  for  SC-cut  crystals  which  have  been  finished  to 

very  loose  manufacturing  tolerances.  The  finished  frequency  and 

the  turning  point  temperature  of  any  ECXO  crystal  can  vary 

considerably  from  the  nominal  design  values  without  any 

degradation  in  accuracy  of  performance  (for  example,  an  ECXO  with 

an  oscillator  frequency  offset  as  large  as  ±  1  x  lO”^  can  be 

.0 

adjusted  to  within  ±  1  x  10  of  the  correct  clock  frequency — a 
compensation  fact<)r  that  is  utterly  unattainable  with  analog 
compensation  methods) . 

With  the  externally  compensated  crystal  oscillator  approach, 
the  oscillator  can  be  expressly  designed  for  maximum  performance 
(highest  Q,  highest  stability,  etc.)  without  any  requirement  for 
the  inclusion  of  a  varactor  or  other  element  that  might  affect 
ultimate  stability. 

The  ECXO  concept  need  not  be  restricted  to  temperature 
compensation  (although  this,  generally  will  be  a  major 
objective)  .  The  compensation  can  be  extended  to  other 
measureable  or  predictable  parameters.  For  example,  the  quartz 
crystal  need  not  be  manufactured  to  a  tight  turnover  temperature 
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or  finishing  tolerance  (an  advantage  if  SC-cut  overtone  crystals 
are  to  be  manufactured  at  a  reasonable  cost)  .  Any  fixed 
frequency  offset#  either  initially  or  observed  at  a  later 
calibration  date#  can  be  readily  removed. 

Finally#  the  long-term  frequency  drift  (aging)  observed  in 
high  stability  quartz  resonators  frequently  proceeds  at  a  nearly 
uniform#  systematic  rate.  ’  This  leads  to  the  possibility  that  a 
microprocessor  ECXO  can  be  programmed  or  calibrated  to  provide 
some  measure  of  compensation  for  systematic  aging. 

The  remarkable  and  continuing  advances  in  solid  state 
technology  and  microcomputer  components  opens  up  the  feasibility 
and  practicality  of  a  moderately  sophisticated  ovenless  ECXO 
system  capable  of  reliable#  accurate  operation  at  very  low  power 
(say#  less  than  20  milliwatt) .  The  power  dissipation  in  CMOS 
integrated  circuits  is  the  sum  of  the  quiescent  (d.c.)  and  the 
dynamic  (a.c.)  components.  The  quiescent  component  becomes 
significant  only  under  conditions  of  high  temperature  and  high 
voltage;  for  example#  the  quiescent  dissipation  of  a  complex  MSI 
device  such  as  the  RCA  C04516B  Presettable  Up/Down  Counter 
approaches  1  milliwatt  at  -t-SS^C  with  a  supply  voltage  of  -t-SVdc. 
The  dynamic  dissipation#  on  the  other  hand#  of  this  counter  has  a 
typical  value  of  roughly  10  milliwatt  at  a  clock  frequency  cf 
SHHz  for  the  same  4-5Vdc  supply. 

The  dynamic  power  in  CMOS  circuits  is  roughly  proportional 
2 

to  the  product  f  x  V  #  where  f  is  the  input  frequency  and  V  is 
the  supply  voltage.  Thus,  a  5:1  reduction  in  power  dissipation 


can  be  achieved  if  a  digital  counter  cc  CMOS  microprocessor  is 
operated  at  1  MHz  instead  of  5  MHz.  Operation  at  reduced  supply 
voltages  also  lowers  power  consumption/  and,  consequently, 
deserves  consideration.  The  higher  power  consumption  will 
undoubtedly  occur  at  the  input  digital  counter/frequency  divider 
stage  in  the  C-mode  channel;  the  B-mode  frequency  channel  can  be 
heterodyned  to  the  lower  difference  frequency  (B-C)  before  any 
digital  frequency  division  process  need  occur,  thereby  lowering 
its  power  requirements.  Once  initial  frequency  division  (say, 
by  a  factor  of  16  or  more)  has  been  achieved,  the  remaining 
digital  circuits,  including  the  microprocessor  and  PROM  memory, 
will  not  add  significantly  to  the  power  budget. 

2.2  Externally  Compensated  Crystal  Oscillator  (ECXO)  Development 
The  initial  Tracor  ECXO  is  shown  in  simplied  block  diagram 
in  Figure  2-1.  This  unit  used  an  SC-cut  crystal  operating  at 
dual  C-mode  (10.23  MHz;  3rd  overtone)  and  B-mode  frequencies. 

The  C-mode  frequency  is  passed  through  a  digital  phase  shifter  (a 
programmable  frequency  divider  which  is  controlled  by  the  MPU 
microprocessor)  and  is  further  divided  to  provide  time-corrected 
1  pps  and  1  kHz  output  signals.  The  temperature  sensitive  B-mode 
frequency,  f^^,  is  down  converted  to  the  difference  frequency 
(f^-f^)  by  the  digital  mixer  measurement  of  > 

in  place  of  can  be  performed  with  a  simpler  lower  frequency 
counter,  thereby  enabling  a  saving  in  power  consumption.  The 
time  base  for  the  frequency  counter  is  derived  by  division  from 
the  time-corrected  C-mode  frequency.  The  measured  value  of 
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is  used  by  the  microprocessor  to  adjust  the  C-mcde 
frequenc^i  in  accordance  with  a  previously  established  calibcaticn 
curve . 

The  time/frequency  adjustment  is  implemented  by  adding  or 
deleting  counts  in  the  digital  phase  shifter.  The  size  of  each 
step  is  roughly  ±2  microseconds  (the  basic  10.23  MHz  oscillator 
frequency  is  divided  down  to  501.2  kHz  because  of  frequency 
limitations  of  the  microprocessor  (MPU)  used  in  the  unit) .  A  ±2 
microsecond  correction  can  be  made  aach  millisecond;  accordingly, 
the  system  is  able  to  handle  frequency  offsets  as  large  as  ±2  x 
10"^. 

It  must  be  noted  that  a  digital  ECXO  of  this  simple  type 
produces  discrete  phase  steps,  rather  than  a  continuous  phase 
rate,  in  the  output  waveforms.  A  ±2  microsecond  quantization 
error  will  be  negligible  in  most  real  time/clock  applications. 

It  will  affect,  however,  the  required  length  of  "calibration” 
time  to  achieve  a  specified  clock  frequency  accuracy  (e.g., 
approximately  100  sec.  of  time  is  required  for  syntonizaticn  tc 
±2  X  10“®,  etc . ) . 

The  Tracer  ECXO  includes  microprocessor  hardware  and 
software  features  to  facilitate  automatic  synchronization  and 
frequency  re-calibration  of  the  system  from  an  external  1  Hz 
reference  signal  pulse.  It  is  merely  necessary  to  connect  a  1 
pps  time  standard  to  the  unit.  The  microprocessor  accepts  an 
incoming  signal  as  a  "valid"  calibration  source  only  if  it  meets 
certain  criteria  with  respect  to  amplitude,  frequency  tolerance, 


and  absence  of  jitter,  thereby  avoiding  any  possibility  that  the 
unit  will  attempt  to  synchronize  on  a  noisy  or  incorrect  signal. 
The  software  is  designed  such  that  the  frequency  calibration 
accuracy  continues  to  be  refined  throughout  the  time  interval 
that  the  calibration  source  is  connected  to  the  unit. 

The  breadboard  ECXO  uses  the  Motorola  6802  CPU 
microprocessor  (with  64  bytes  of  internal  RAM) ,  in  combinaticn 
with  a  Motorola  6840  timer  chip,  for  the  major  control,  logic, 
and  counter/timing  functions.  The  instruction  program  and  the 
compensation  look-up  table  values  are  stored  in  a  standard  EPRCM 
memory  (a  2716  NHOS  2048  x  8-bit  OV  eraseable  PROM) .  Several  low 
power  Schottky  integrated  circuits  (e.g.,  74LS74,  74LS390,  etc.) 
are  also  used  for  auxiliary  gating  and  buffer  amplifier  purposes. 

The  required  temperature  compensation  pattern  for  each 
SC-cut  crystal  was  individually  determined  by  means  of 
calibration  in  a  temperature  chamber.  To  collect  the  raw  data 
from  which  the  look-up  table  was  prepared,  simultaneous  readings 
of  C-mode  and  B-mode  (or,  alternatively,  frequencies  were 

taken  as  the  temperature  of  the  environmental  test  chamber  was 
cycled  over  the  full  temperature  range  (the  exact  temperature 
need  not  be  known  since  temperature  only  enters  parametrically 
into  the  frequency  ratio)  .  For  the  case  of  the  two  ECXC 

units,  it  was  found  that  the  raw  data  for  each  SC  crystal 
approximated  a  cubic  equation.  A  general  purpose  mini-computer 
was  used  for  verifying  this  and  for  determining  the  coefficients 
of  the  resultant  cubic  equation.  The  individual  compensation 


values  to  be  programmed  into  the  PROM  look-up  table  were  then 
determined,  again  with  a  general  purpose  computer,  so  as  to 
properly  match  the  cubic  equation  at  selected  intervals  along  the 
temperature  range. 

The  com.pensation  algorithm  within  ECXO  memory  provides  a 
linear  interpolation  procedure  between  the  discrete  points  stored 
in  the  table.  The  table  incorporates  some  200  points,  each  16 
bits  in  length. 

Computer  simulation  of  the  compensation  algorithm  was  used 

to  demonstrate  that  the  total  compensation  program  would  operate 

satifactorily,  without  introducing  any  significant  systematic  or 

random  error  in  either  clock  frequency  or  timing  over  extended 

intervals.  These  computer  simulation  runs  demonstrated  that  the 

compensation  algorithm  in  itself  would  introduce  considerably 

less  than  1  millisecond  time  error  over  a  24-hour  operating 

interval  in  a  typical  temperature  environment. 

Initial  testing  of  the  two  demo  units  under  fluctuating  room 

ambient  conditions  indicated  that  a  timing  accuracy  of  better 

than  ±1  millisecond  could  generally  be  maintained  for  1-3  days 

after  time  sychronization  and  frequency  re-calibration. 

Re-calibration  was  performed  by  means  of  a  laboratory  standard  1 

pps  signal  having  an  accuracy  of  ±1  x  lO”^*^  or  better.  In  sore 

q 

instances  a  systematic  drift  of  several  parts  in  10'  was  observec 

throughout  the  test  period;  in  other  runs  the  time  error  appeared 

to  follow  a  random  walk  pattern,  with  the  apparent  frequency 

—  8 

fluctuating  over  a  range  of  approximately  ±2  x  10  during  a 


several-hour  period. 

Subsequent  environmental  chamber  testing  of  one  of  the  units 
several  months  later  showed,  however,  an  inexplicable  frequency 
offset  of  perhaps  1  x  lO”^.  The  question  then  arose  whether  the 
deviation  was  due  to  an  intrinsic  change  in  the  B-mode  or  C-mode 
resonant  frequency  of  the  quartz  crystal  or  to  a  change  in  sore 
circuit  component.  To  resolve  this  crucial  question  it  was 
decided  that  basic  measurements  of  the  B-mode  and  C-mode 
frequency  characteristics  of  SC-cut  resonators  were  prerequisite 
to  the  ECXO  program. 


SECTION  3 


INVESTIGATION  OF  DUAL  MODE  (B  &  C)  FREQUENCY 
CHARACTERISTICS  OF  SC-CUT  CRYSTALS 

3.1  INTRODUCTION 

The  early  ECXO  used  a  modification  of  the  familiar  Butler 
oscillator  circuit  for  simultaneous  excitation  of  the  B-mode  and 
C-mode  frequencies.  There  is  an  inherent  problem  in  any  active 
feedback  type  of  oscillator  that  attempts  to  generate  two 
independent  frequencies  from  a  dual  mode  crystal.  Frequency 
selective  networks  must  be  employed  to  separate  the  modes  and  to 
prevent  undesirable  interaction  between  the  individual  frequency 
modes.  The  electronic  parameters  of  the  components  used  in  the 
frequency-selective  feedback  loops  will  show  some  residual  drift 
with  time  and  changes  in  environmental  conditions.  Consequently 
the  long-term  stability  and  aging  characteristics  of  any 
conventional  oscillator  used  in  a  dual  B-raode  and  C-mcde 
resonator  application  can  be  questioned. 

3.2  THE  PHASE  MODULATED  FREQUENCY  DOMAIN  REFLECTOMETER 

A  review  was  therefore  made  of  alternative  methods  for 
laboratory  investigation  of  the  B-raode  and  C-mode  frequencies  of 
precision  SC-cut  resonators.  As  a  goal,  such  instrumentation 
should  be  able  to  provide  a  significantly  higher  level  of 
measurement  precision  (say,  by  a  factor  10  or  more)  than  that 
required  in  an  operational  ECXO.  An  instrumentation  system 
employing  a  unique  phase  modulated,  frequency-lock  r ef lectcm.eter 
was  selected;  Figure  3-1  shows  the  basic  concept.  The  crystal  i 
used  here  as  a  passive  resonator  in  a  frequency-lock  servo 


oscillator  system.  The  advantages  of  this  systeri.  include:  (1) 
the  crystal  resonator  may  be  remotely  located  (e.g.,  in  an 
environmental  test  chamber)  away  from  the  rest  of  the 
instrumentation;  (2)  significantly  improved  immunity  to  component 
parameter  values;  (3)  the  potential  for  almost  complete  freedom 
from  interaction  between  the  C-mode  and  B-mode  portions  of  the 
ECXO;  (4)  the  ability  to  achieve  accurate,  repeatable  crystal 
measurements  and  a  quality  of  temperature  compensation  that  is 
limited  only  by  uncertainties  intrinsic  to  the  resonator;  and  (5) 
the  ability  to  maintain  a  specified,  constant  level  of 
independent  excitation  for  C-mode  and  B-mode  frequencies  without 
a  requirement  for  AGC. 

An  external,  independent  voltage-controlled  oscillator 
generates  a  carrier  frequency,  f  ,  in  the  vicinity  of  the  crystal 
resonant  frequency,  A  small  phase  modulation,! a  0  ,  is  applied, 
and  the  resultant  phase  modulated  signal  is  used  to  interrogate 
the  crystal  resonator.  The  reflected  signal  from  the  resonator, 
after  passing  through  an  isolation  amplifier,  is  amplitude 
detected  by  a  diode  and  capacitor  filter.  An  amplitude 
modulation,  at  the  modulation  oscillator  frequency,  is  observed 
in  the  detector  output  whenever  the  carrier  frequency  differs 
from  the  resonant  frequency,  f^,  of  the  quartz  resonator.  The 
sense  and  magnitude  of  this  detected  amplitude  modulation  depends 
tipon  the  frequency  deviation  of  the  voltage-controlled  oscillator 
from  the  center  of  crystal  resonance. 

A  phase  sensitive  detector  is  used  to  provide  a  d-c  error 


signal  proportional  to  the  frequency  offset  of  the  oscillator. 
This  error  signal;  after  suitable  low  pass  filtering,  is  used  to 
servo  control  the  oscillator,  thereby  frequency-locking  the 
oscillator  to  the  crystal  resonant  frequency. 

It  will  noted  that  Figure  3-1  utilizes  a  directional  coupler 
(or  four-port  hybrid  transformer)  for  coupling  the  phase 
modulated  input  signal  to  the  quartz  resonator.  The  use  of  the 
hybrid  transformer  enables  good  isolation  to  be  achieved  betv/een 
the  input  to  the  crystal  and  its  reflected  output  signal. 

Dual  mode  operation  from  a  single  resonator  can  be  obtained 
by  utilizing  separate  frequency-lock  oscillators  for  the 
individual  C-mode  and  B-mode  frequencies  (Figure  3-2) . 

Interaction  between  the  two  modes  can  be  eliminated  by  using 
separate,  non-r elated  modulation  frequencies  for  the  phase 
modulation/demodulation  processing,  together  with  good  low  pass 
filtering  in  the  servo  control  loop. 

A  better  understanding  of  the  operating  principles  of  the 
modulated  system  can  perhaps  be  otained  from  the  vector  diagrams 
of  Figure  3-3.  and  R2  represent  the  successive  vectors 
associated  with  the  incident  signal  as  the  result  of  the  square 
wave  phase  modulation  process;  S^,  and  C  represent  the 
equivalent  sidebands  and  carrier  components  for  this  phase 
modulated  signal.  The  corresponding  reflected  components  are 
shown  on  the  right.  The  crystal  resonator  unit  appears  as  a  veri- 
high  impedance  (highly  reactive)  for  all  frequency  components 
even  slightly  removed  from  resonance;  accordingly,  the  reflected 


sideband  cor.ponents,  and  S'2»  are  virtually  equal  to  the 
incident  components,  and  S21  and  maintain  the  same  180° 
relationship  to  each  other.  It  is  the  carrier  com.ponent  that 
changes  in  magnitude  and  phase  as  the  crystal  resonant  frequency 
is  approached.  At  resonance,  there  is  likely  to  be  a  residual 
reflected  component,  C  (as  the  result  of  any  nominal  mismatch 
between  the  ESR  of  the  crystal  and  the  hybrid  termination 
resistance  •  However,  this  C  component  (at  resonance)  is 
orthogonal  to  the  sideband  component  vectors  so  that  the 
resultant  vectors,  P.'^  and  R'2»  are  equal  in  magnitude.  Away 
from  resonance,  the  reflected  carrier  component  (C  '  or  C'') 
increases  in  magnitude  and  undergoes  a  large  phase  shift; 
accordingly,  the  two  vectors,  R’^^  and  R'2,  are  no  longer  equal  in 
magnitude  and,  consequently,  an  error  signal  is  generated  in  the 
synchronous  demodulation  process. 

The  length  of  the  interconnecting  cable  between  the  quartz 
resonator  and  the  balance  of  the  instrumentation  does  not 
materially  influence  the  measured  resonance  frequency  if  the 
hybrid  bridge  is  properly  terminated  to  match  the  cable 
impedance.  An  ideal,  non-dispersive  cable  would  affect  the 
carrier  and  sideband  frequency  components  in  an  identical  manner, 
thereby  producing  no  change  in  the  relative  magnitude  or  phase  of 
the  detected  vector  components.  Furthermore,  for  moderate 
lengths  of  commonly  available  transmission  cables,  the  dispersion 
effect  is  quite  small  and  can  generally  be  neglected.  This  has 
been  confirmed  in  simple  tests  wherein  an  additional  5C  feet  of 


RG-58U  cable  has  been  inserted  between  the  resonator  and  the 

ref lectoip.eter  bridge;  without  using  extreme  care  in  adjusting  the 

-c 

hybrid  termination,  a  frequency  deviation  of  less  than  5  x  1C 
was  observed. 

3.3  DATA  ACQUISITION  AND  ANALYSIS  INSTRUMENTATION 

Figure  3-4  shows  the  measurement  system  used  during  the 

early  phases  of  the  project.  A  brassboard  dual  channel 

ref lectometer  enabled  the  measurement  of  the  B-mode  and  C-moce 

resonant  frequencies  of  the  SC-cut  quartz  crystal  to  a  precision 
.a 

of  ±2  X  10  .  The  crystal  resonator,  together  with  the  hybrid 

transformer,  was  remotely  located  within  the  environmental  test 
chamber.  Standard  RG-58U  coaxial  cable,  approximately  5  meters 
in  length,  linked  the  crystal  resonator  to  the  ref lectometer . 

A  timer,  with  a  total  cycle  time  adjustment  between  1-10 
minutes,  controlled  the  multiplexer  (MPX) .  During  each  cycle  the 
two  R-F  frequencies  were  sequentially  switched  in  the  pattern:  E 
(for  1  second  duration);  C  (for  3  seconds);  and  B  again  (1 
second) .  The  first  and  second  B-frequency  values  were  averaged 
together  to  provide  an  "interpolated"  B-mode  reading  matching  the 
time  frame  of  the  intermediate  C-reading.  Simultaneous  D  and 
C-mode  readings  with  two  counters  would  have  obviated  this 
requirement  for  a  multiplexer  and  dual  B-mode  readings. 

The  HP-5335A  is  a  highly  versatile  counter  that  provides  a 

-9 

precision  and  resolution  of  ±  2  x  10  fo-r  a  gating  interval  of  1 
second.  The  3-second  interval  for  the  more  critical  C-mcoe 


FIGURE  3-4 


frequency  enables  a  counter  measurement  precision  in  excess  of  1 

-9 

X  10  .  All  data  throughout  the  project  were  taken  relative  to  a 

laboratory  reference  standard  having  an  accuracy  of  +2  x  IG  .  or 
better  (using  periodic  calibrations  based  on  Loran  C  signal 
reception) . 

In  the  early  work  the  individual  frequency  readings  were 
simply  outputted  from  the  counter  onto  a  paper  tape  printer. 

After  completion  of  a  full  temperature  run^  the  tape  data  were 
visually  read  and  manually  keyed  into  the  HP-87  computer  for 
permanent  storage  and  analysis. 

Discrete  frequency-temperature  data  of  this  type  can  be 
analyzed  in  various  ways.  The  approach  which  appears  to  be  most 
useful  involves  determination  of  the  frequency  deviation  of  the 
data  from  an  assumed  power  series  polynomial  of  the  Mth  degree. 

It  is  known  that  the  frequency-temperature  dependency  of  the 
C-mode  frequency  can  be  reasonably  approximated  with  a  cubic 
equation.  However f  with  the  powerful  computer  techniques 
available  today,  it  is  readily  possible  to  utilize  a  higher 
degree  polynomial  expansion,  if  this  will  provide  more  useful 
interpretation  of  the  data. 

All  of  the  data  in  the  subsequent  section  are  presented  ir. 
terms  of  the  residual  frequency  deviation  of  the  individual  data 
points  from  an  "ideal"  or  predicted  power  series  polynomial  of 
the  form 

f  -  *0  *  V  ♦  * . *  v" 

where  P  is  the  predicted  frequency  and  x  is  a 


temperature-dependent  frequency  variable  (e.g.,  the  B-C  frequency 
difference)  that  forms  the  basis  for  the  temperature  compensaticn 
technique.  The  individual  polynomial  coefficients  in  the  series 
can  be  computed  by  employing  standard  regression  techniques  on  a 
selected  set  of  C-mode  and  (B-C)  frequency  data  pairs.  (The 
software  program  "REGR/2107'  in  Appendix  B  of  this  Final  Report 
utilizes  this  capability) . 

It  might  be  noted  that  deriving  a  high  order  polynomial 
(e.g.f  N  >  5  or  higher)  from  a  large  number  of  data  pairs 
requires  considerable  computer  time  and,  in  addition,  places  some 
burden  on  computer  accuracy  if  the  magnitude  of  the  variable,  x, 
is  too  large.  Consequently,  the  B-C  frequency  value  is  generally 
offset  and  reduced  by  an  arbitrary  constant  in  establishing  the 
variable  x  used  in  Equation  1. 

Figure  3-5  shows  the  improved  automatic  data  logging  system 
employed  during  the  later  stages  of  the  project.  This 
configuration  utilizes  the  HP-IB  (IEEE-488)  Interface  Bus  for 
system  control  and  transfer  of  data.  The  HP-87  computer  forms 
the  "controller"  for  the  system  and,  in  addition,  can  be  used  for 
preliminary  plotting  and  analysis  of  the  raw  data  on  a  real-time 
basis.  For  this  purpose  it  is  merely  necessary  that  the  operator 
insert  an  appropriate  set  of  polynomial  coefficients  for  the 
particular  crystal  (s)  under  test.  The  computer  will  then  print 
out  and  graphically  plot,  on  real-time  basis  throughout  each 
temperature  run,  the  frequency  deviation  of  the  successive  C-mode 
frequency  values  from  the  predicted  polynomial.  The  program 


listing,  "DATA/PLOT*,  in  Appendix  C  gives  the  actual  program  used 
for  (a)  controlling  the  ref lectoir.eter  and  counter  gating;  (b) 
storing  the  raw  data  sets  on  disk;  (c)  operator  inputting  of  the 
polynomial  coeif f icients;  (d)  hard  copy  printout  of  the 
individual  "raw"  frequency  data  and  time;  (e)  printout  of  the 
frequency  deviation,  L  f/f,  of  the  C-mode  frequency  from  its 
corresponding  predicted  value  based  on  the  assumed  polynomial 
approximation;  (f)  a  plot  of  these  same  frequency  deviation 
values  on  the  CRT  terminal;  and  (g)  on  command,  a  "dump"  of  the 
resultant  CRT  graphics  display  onto  the  printer/plotter. 

3 . 4  TEST  RESULTS 

The  present  investigation  into  the  behavior  and 
interrelationship  of  the  dual  mode  frequency  characteristics  of 
SC>cut  quarts  resonators  has  demonstrated  the  tremendous 
potential  of  the  dual  mode  ECXO  in  low  power,  ovenless  clock 
driver  and  timing  applications.  All  results  obtained  to  date, 
although  limited  to  a  relatively  small  sample  of  SC-cut  cristals 
from  several  suppliers,  indicate  an  attainable  frequency  accuracy 

O 

of  several  parts  in  10  over  the  full  military  temperature  range 
of  -55°C  to  -J-SS^C  (corresponding  to  an  equivalent  clock  accuracy 
of  several  milliseconds/day  or  better).  This  approaches  a 
100-fold  improvement  over  that  obtainable  in  practice  with  a 
varactor-type  TCXO  over  the  sarae  extreme  temperature  range.  It 
must  be  noted,  of  course,  that  the  any  accuracy  value  given  here 
is  based  solely  on  the  intrinsic  stability  of  the  dual  mode 
quartz  resonator  that  forms  the  nucleus  of  the  ECXO;  however, 


with  further  research,  it  should  be  possible  to  develop  suitable 
oscillator  and  associated  microprocessor  circuits  that  do  not 
materially  degrade  the  intrinsic  resonator  performance. 

A  total  of  several  hundred  thousand  data  points,  covering 
more  than  fifty  separate  environmental  test  chamber  runs,  were 
probably  logged  during  this  investigation.  Many  runs, 
particularly  during  the  early  months  of  the  work  program,  were 
duplicated  to  ensure  that  the  instrumentation  was  working 
properly  and  that  valid  results  were  being  obtained. 

Accordingly,  only  a  small,  representative  sample  of  the  more 
recent  test  data,  together  with  a  review  of  the  significant 
conclusions  that  can  be  drawn  from  the  data,  will  be  represented 
here. 

Table  3-1  illustrates  the  format  used  for  print-out  of  the 
raw  data,  together  with  interim  processing  of  these  data,  on  a 
real-time  basis.  The  first  40  lines  of  a  temperature  run  (with 
the  full  run  covering  a  60  hour  cycle)  are  reproduced  here.  Eac 
set  of  data  logging  includes  two  sequential  pairs  of  dual-mode 
frequencies  in  the  sequence  E^,  B^',  C^'.  The  values  ir.  the 

column  TMODEl  are  the  average  of  the  and  readings;  the 
column  TMODE  shows  the  differential  between  these  two  readings 
The  column  CMODE  gives  the  CMODE  readings;  tlie  column  CKCDE 
shows  the  difference  between  successive  rows  of  CMODE  values 
(here  separated  by  a  sampling  time  interval  of  approximately  51 
seconds,  as  evident  by  inspection  of  the  TIKE  column) .  TAPPROX 
approximates  the  apparent  temperature  of  the  quartz  resonator 
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NO. 

TFIOOBl 

-TWOOE 

cnooE 

-CnODE 

TARRROX 

TINE 

-R/R  DEV 

X  KZ 

1 

S47E134.09 

-.37 

4999972.131  . 

♦24.01 

394*2 

t 

S47E13S.39 

-.*8 

4999972. 139 

-.012 

•24.00 

39912 

♦  1.S2E-OS 

♦8.  1*4 

S47E13A.99 

-.80 

4999972.  124 

-.013 

♦23.99 

399*3 

♦1.S1E-08 

-S.  1*5 

4 

S47E13E.E3 

-.89 

4999972. 108 

-.01* 

♦23.98 

39*14 

♦1.84E-08 

-8. 1*7 

3 

S47E140.92 

-1.03 

4999972.089 

-.019 

♦23.9* 

39**3 

♦1.54E-08 

♦9. 1*9 

* 

S47E145.3S 

-1. 19 

4999972.  0*7 

-.022 

♦23.94 

3971* 

♦1.5SE-08 

♦8. 172 

7 

S47E143.9S 

-1.21 

4999972.  043 

-.024 

♦23.92 

397*8 

♦1.55E-08 

♦8. 175 

E 

347E14E.71 

-1.24 

4999972.019 

-.024 

♦23.90 

39818 

♦1.57E-08 

♦9. 177 

* 

S47E191.SS 

-1.47 

4999971.992 

-.027 

♦23. 88 

398*8 

♦1.S6S-08 

♦8. 180 

10 

S47E1SS.3S 

-1.99 

4999971.938 

-.034 

♦23. 89 

39919 

♦1.60E-08 

-8. 184 

11 

347ElS9.Sfl 

-2.21 

4999971.913 

-.043 

♦23.82 

39970 

♦1.61E-08 

♦8. 189 

13 

S47El*4,95 

-2.29 

4999971.8*8 

-.047 

♦23. 78 

40021 

♦1.60E-08 

-8. 194 

13 

3479170.0* 

-2.23 

4999971.822 

-.04* 

♦23.75 

*0073 

♦1.61E-08 

♦9. 199 

14 

S47E173.3* 

-2.  *4 

4999971.772 

-.090 

♦23.71 

40123 

•1.62E-08 

♦9.205 

IS 

S47E1E1.49 

-2.88 

4999971.718 

-.03* 

♦23.6* 

4017* 

♦1.69E-08 

♦8.211 

1* 

S47E1SE.07 

-2.  97 

4999971. *37 

-.0*1 

♦23.61 

40229 

♦1.70E-O8 

♦8.218 

17 

S47E194.74 

-2.94 

4999971.396 

-.0*1 

♦23.3* 

4027* 

♦1.72E-08 

♦8.225 

IE 

S47E301.41 

-3.  18 

4999971.532 

-.0*4 

♦23.31 

4032*  • 

♦1. 7IB-08 

•8.231 

19 

347E20E.aS 

-3.58 

4999971.4*3 

-.0*9 

♦23.4* 

40377 

•1.74E-08 

•9.239 

30 

547E317.09 

-3.78 

4999971.38* 

-.077 

♦23.40 

4042* 

♦1. 79S-Oe 

•9.248 

Cl 

S47S32S. 38 

-3.74 

4999971.303 

-.083 

♦23.3* 

40479 

♦1.74C-08 

•8.256 

5478253.99 

-3.73 

4999971.226 

-.077 

♦23.28 

40530 

♦i.80E-oe 

-9.265 

3478242. 53 

-3.99 

4999971. 143 

-.081 

♦23. 21 

40581 

-1. 845-06 

-9.273 

24 

3478331. *4 

-4.  12 

4999971.058 

-.087 

♦23. 13 

40*32 

♦1.87E-08 

•8.293 

34782*0. 93 

-4.  13 

4999970. 9*8 

-.090 

♦23. 08 

40*83 

♦1 • 87S— 08 

♦9. 292 

3* 

3479270. 17 

-3.99 

4999970.881 

-.087 

♦23.01 

40733 

•1.91E-08 

-9.301 

37 

5479279.28 

-4.37 

4999970. 790 

-.091 

♦22.99 

40785 

♦1. 9IE-08 

•a. 311 

3E 

5479289.40 

-4.  72 

4999970.  690 

-.  100 

♦22.87 

4083* 

- 1 . 94E-08 

-9.321 

29 

3478300. 09 

-4.74 

4999970.583 

-.  103 

♦22.79 

4088* 

♦ 1 . 9*8-08 

-9.332 

30 

3478310. 70 

-4.  *2 

4999970. 482 

-.  103 

♦22.72 

40937 

•2. OOE-OB 

— S . —4^ 

31 

3478321.01 

-4.57 

4999970.381 

-.  101 

♦22. *4 

40988 

♦2.05E-O8 

-8.333 

5478331.4* 

-4.  74 

4999970. 273 

-.  10* 

♦22.97 

41040 

-2. 05E-03 

•a.  364 

5478342. 15 

-4.74 

4999970. 1*8 

-.  107 

♦22.49 

41090 

♦2.08E-O8 

-a. 374 

34 

5478332. 76 

-4.60 

4999970. 0*1 

-.  107 

♦22.41 

41141 

-2. lOE-09 

-9.335 

39 

54783*2. 99 

-4.38 

49999*9. 9*0 

-.  101 

♦22.34 

41192 

♦2. laE-Oe 

-a.  395 

3* 

3478373. 02 

-4.  75 

49999*9. 3»3 

-,  107 

•22.27 

41243 

•2.  13E-oa 

-a. 406 

37 

3478383.93 

4.98 

49999*9. 742 

-.111 

♦22. 19 

4  1294 

♦2. 20E-ce 

-9.41' 

3B 

5478395. 14 

-4.91 

49999*9. *25 

-.117 

♦22. 10 

41345 

-2.  20E-08 

-8. *28 

39 

347840*. 1 1 

-4.73 

49999*9.313 

-.  1 12 

♦22.03 

4139* 

•2. rsE-oa 

-a. 439 

40 

947841*. *3 

-4.  73 

49999*9. 401 

-.  112 

♦21.95 

4144* 

•2. I4E-08 

♦8.  450 

TYPICAL  PRINTOUT  OF  B  &  C-MODE  DATA 
ON  REAL-TIME  BASIS 
TABLE  3-1 
3-1^ 


under  test  (based  upon  an  assumed  linear  relationship  between 
B-mode  frequency  and  temperature  of  -138Hz/°C  over  the  full 
operating  range) . 

The  tabulated  frequency  deviation  values,  i  F/F,  are  shown 
relative  to  a  5th  degree  polynomial.  The  polynomial 
coefficients,  once  derived  from  an  initial  "calibration"  run,  ca 
be  utilized  in  all  subsequent  runs  to  measure  hysteresis,  aging 
effects,  etc.  The  last  column  shows  the  offset  B-C  frequency 
value  (in  kHz)  that  is  used  as  the  independent  variable  in  the 
polynomial  frequency  correction  expression. 

Certain  initial  information  must  be  provided  by  the 
operation  prior  to  a  new  run  (e.g.,  type  of  crystal,  incident 
crystal  power,  planned  temperature  profile,  etc.).  This 
ancillary  information,  as  shown,  is  also  stored  on  each  disk  and 
is  used  to  aid  in  subsequent  identification  and  analysis. 

In  addition  to  the  hard-copy  digital  printout,  the  computer 
terminal  provides  a  real-time,  ongoing  display  of  the  frequencv 
deviation  of  the  individual  readings  as  they  are  accumulated. 
This  visual  display  can  then  be  "dumped"  onto  a  printer/plotter 
upon  command.  Such  a  plot  is  shown  in  Figure  3-6.  This 
represents  a  total  of  4,000  readings,  taken  over  a  60  hour 
period,  at  oven  intervals  of  less  than  0.1°C.  The  lighter,  uppe 
trace  indicates  data  taken  during  the  decreasing  temperature 
portion  of  the  cycle;  the  heavier  trace  (formed  by  the  broad 
sympbol  'o',  rather  than  a  simple  indicates  increasing 

crystal  temperature. 


PLOT  OF  FREQUENCY  DEVIATION  FOR  SC-CUT  CRYSTAL 
(COLORADO  CRYSTAL  CO.  S/N  1)  DURING  SLOW  TERMPERATURE 
CYCLING  (l”C/12  MINUTE).  UPPER  TRACE:  DATA  DURING 
DECREASING  TEMPERATURE;  LOWER  TRACE:  INCPvEASING 
TEMIU'PvATURE 


Several  significant  features  can  be  noted  in  these  plots. 

Even  at  the  relatively  low  temperature  rate  of  ±5°/hour/  there  is 

a  clear-cut  displacement  between  the  increasing  and  decreasing 

temperature  curves:  roughly,  a  frequency  differential  of  ±1.5  x 
—  8 

10  is  seen  between  the  negative  rate  and  positive  rate 
sections.  There  is  also  an  observable  activity  dip  in  the 
B-mode;  for  this  particular  crystal,  the  transition  occurs  at 
approximately  +68°C.  Similar  B-mode  activity  dips  have  been 
observed  in  the  other  crystals  that  have  been  tested.  In  some 
cases,  it  has  not  been  possible  to  maintain  satisfactory 
operation  over  the  entire  -55°C  to  +85°C  range;  however,  for  this 
crystal  (Colorado  Crystal  Co;  S/N  1)  the  dip  shows  up  simply  as 
an  abrupt  transition,  with  a  frequency  deviation  step  of  less 
than  3  X  10  in  the  apparent  C-mode  compensation  term. 

Figures  3-7  and  3-8  show  comparable  performance  at  higher 
oven  rates.  Figure  3-7  si  jws  one  complete  temperature  cycle  at 
an  oven  rate  of  ±l'^C/3  minute  (roughly,  quadruple  the  rate  used 
for  the  Figure  3-6  run).  However,  there  is  only  a  moderate 
increase  in  the  separation  of  the  upward  and  downward  temperature 
regions . 

Figure  3-8  shows  the  cluttered  "real  time"  plot  that  results 
when  all  data  from  a  complex  temperature  run  are  displayed  or.  a 
single  plot.  In  this  instance,  the  quartz  resonator  was  cycled 
over  two  separate  low  temperature  and  high  temperature 
sub-regions  as  well  as  an  additional  two  cycles  over  the  full 
-55°C  to  +95°C  range.  The  computer  program  is  set  up  sc  that,  in 
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subsequent  analysis,  selected  segments  of  the  stored  raw  data  can 
be  retrieved  and  analyzed.  This  capability  is  shown  in  Figure 
3-9.  The  left-hand  hysteresis  loop  shows  an  intermediate  segment 
in  which  the  oven  temperature  vas  increased  and  then  decreased  so 
that  the  crystal,  while  at  a  temperature  of  -32°C  from  a  previous 
temperature  cycle,  was  raised  to  roughly  -5°C  and  then  returned 
to  -32°C  (over  a  time  interval  of  1  hour).  The  right-hand 
pattern  shows  similar  resonator  behavior  in  the  high  temperature 
region  between  +60°C  and  +8S°C:  the  upper  dotted  trace  shows  the 
frequency  deviation  during  the  temperature  decrease,  and  the 
points  with  the  square  symbol  show  data  taken  during  increasing 
temperatures.  The  step  transition  due  to  the  activity  dip  is 
clearly  evident  in  each  trace. 

The  most  striking  feature  of  the  dual  mode  frequency 

behavior  is  the  high  stability  and  repeatability  of  the  frequency 

deviation  patern  for  successive  runs  under  similar  test 

conditions.  There  is  also  clear  evidence  of  a  repeatable 

hysteresis  or  thermal  shock  effect  during  temperature  cycling. 

The  test  results  indicate  that  simple  temperature  compensation 

(e.g.,  without  additional  compensation  for  the  rate  of 

temperature  change)  cannot  yield  the  highest  accuracy.  For 

—8 

example,  to  achieve  1  x  10  accuracy  it  will  be  necessary  that 
some  form  of  thermal  rate  compensation  be  employed.  If  the 
functional  relationship  can  be  determined,  either  or.  an  empirical 
or  physical  basis,  any  required  additional  computational  task  can 
be  readily  accomplished  by  simple  microprocessor  means. 


For  static  temperatare  compensation,  the  computed  frequency 
compensation  term,  P(TEMP),  a  function  of  temperature,  is  assumed 
to  approximate  a  polynomial  power  expression  of  the  form 

P(TEMP)  ■  Aq  +  Aj^x  +  AjX^  + . +  AjjX^.  (Eq.  lA) 

Thermal  transient  compensation  can  be  included  by  modifying  this 
expression  to  form  a  new  compensation  factor, 

P’(TEMP)  ■  P(TEMP)  +  K  II  K(2)*(T2  -  25)1 

where  K  and  R(2)  are  thermal  rate  coefficients  and  (T2  -  T^) 
represents  the  change  in  crystal  temperature  over  the  time 
interval  (t2  -  tj^) .  The  K  factor  introduces  a  frequency 
compensation  value  that  is  directly  proportional  to  the  thermal 
rate  at  any  point  in  the  temperature  cycle;  the  R(2)  factor,  on 
the  other  hand,  introduces  a  secondary  correction  that  also 
controls  the  rate  compensation  as  a  function  of  crystal 
temperature  (thereby  permitting  differential  adjustment  of  the 
rate  compensation  over  the  low  temperature  and  high  temperature 
regions) . 

A  limited  number  of  data  analysis  runs  were  made  to 
determine  a  pair  of  suitable  values  for  both  the  R  and  the  R(2) 
coefficients.  For  example.  Figure  3-10  shows  the  frequency- 
deviation  relative  to  a  derived  fifth  order  frequency-temperature 
polynomial  for  a  dual  mode  SC-cut  crystal  without  rate 
compensation.  The  upper  trace  (small  dots)  represents  data  taken 
during  decreasing  temperatures;  the  lower  trace  (sm.all  open 
squares)  indicates  corresponding  points  for  increasing 
temperatures.  The  environmental  chamber  was  programmed  to  give  a 


‘^2  -  ^1’ 

(t2  -  t^) 


(Eq.  2) 


3-20 


ttt  STuNOftRO  DEVIATION  ('•F/F)  «  1 . 94S70276SS1E-B 

tf«  AVE  ABSOLUTE  DEVIATION  (F/F)-  1 . S021 1 70S705E-8 

ttt  AVE  FREQ  DEVIATIQNC  F/F'  1 . 36177697701E-B 

ttAuCUMULATED  tine  RATE  ERROR ( M I LL I SECONDS / DA Y > •  1 . 17674BI0814 

thermal  shock  coefficient (F/F  per  deg  C/ SECOND) I  0 

K(2>«  0 


temperature  ramp  at  an  approximate  rate  of  ±0.5°C  per  minute.  A 
frequency  deviation  differential  in  excess  of  2  x  10  is  obvious 
between  the  decreasing  and  increasing  portions  of  the  temperature 
cycle.  Figure  3-11  shows  the  same  data  but  with  a  thermal  rate 
compensation  factor  added  to  the  static  polynomial  temperature 
compensation.  It  will  be  noted  that  the  previous  differential 
offset  has  been  markedly  reduced. 

Figure  3-12  and  Figure  3-13  show  similar  results  for  another 
run>  several  days  later,  at  a  higher  thermal  rate  (slightly  in 
excess  of  ±l®C/minute) .  Figure  3-12  shows  test  results  with  only 
static  temperature  compensation  employed.  Figure  3-13,  with  a 
noticeably  smaller  deviation  differential,  is  based  on  the  same 
rate  compensation  factor  used  in  Figure  3-11  (namely,  K  »  -1  x 
lO”^  per  ®C/second,  and  K(2)  «  -0.005).  Inspection  of  Figure 
3-13  indicates  a  small  measure  of  over compensation  in  the  thermal 
rate  factor.  A  value  in  the  proximity  of  K  ■  -0.9  x  10“®  would 
perhaps  have  provided  better  rate  compensation  over  the  full 
temperature  cycle  of  this  particular  run.  Figures  3-14  and  3-15, 
for  the  case  of  an  intermediate  thermal  cycle  over  the  low 
temperature  region  of  -8°C  to  -32®C  and  return,  again  illustrate 
the  advantages  of  incorporating  thermal  rate  compensation.  V^e 
again  see  evidence  of  overcompensation  if  we  employ  the 
"baseline"  values  of  K  ■  -1  x  lO”^  and  K(2)  ■  -0.005  that  had 
been  derived  from  the  data  in  Figure  3-11. 

From  these  data  it  can  be  concluded  that  the  linear  form  of 
thermal  rate  compensation  employed  in  these  tests,  although 
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offering  a  very  useful  and  significant  improvement  in  accuracy, 
does  not  provide  complete  compensation  over  varying  envirorniental 
conditions.  A  non-linear  expression  (e.g.,  one  using  the  square 
root  of  the  measured  thermal  rate)  or,  alternatively,  inclusion 
of  a  thermal  acceleration  term  might,  on  an  empirical  basis,  give 
better  results  over  a  wider  range  of  environmental  conditions. 

A  portion  of  the  observed  frequency  displacement  during 
temperature  cycling  may  be  due  to  hysteresis;  however,  it  is 
difficult  to  separate  hysteresis  effects  from  the  residual 
frequency  differencies  that  are  observed  at  the  lowest  thermal 
rates  (±5°C  per  hour)  used  in  some  of  the  environmental  testing. 

Figure  3-16  illustrates  the  comparative  frequency  deviation 
pattern  of  an  SC-cut  lOMHz  crystal  cut  to  have  a  lower  turnover 
temperature  in  the  vicinity  of  +5°C  for  the  C-mode  {3rd  overtone) 
frequency.  The  dotted  trace  shows  the  response  for  decreasing 
temperatures;  the  heavier  trace  shows  the  increasing  temperature 
portions  of  the  thermal  ramp  cycle  (at  the  slow  rate  of 
±0 .25°C/minute) .  This  pattern  differs  from  the  previous  plots  in 
that  the  thermal  transient  behavior  appears  to  reverse  polarity 
in  the  neighborhood  of  the  turning  point:  at  elevated 
temperatures,  the  observed  deviation  for  decreasing  temperatures 
falls  below  the  readings  for  increasing  temperatures.  Also,  it 
appears  that  the  thermal  shock  (or  hysteresis)  effect  is 
relatively  small  throughout  the  low  temperature  region. 
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It  should  be  recognized  that  plotting  of  the  data  in  terr.s 
of  the  frequency  deviation  from  any  particular  polynomial  or 
other  analytic  "baseline  reference"  is  an  expedient  that 
facilitates  data  analysis  and  interpretation.  The  observed 
deviation  and  structure  of  the  "frequency  deviation"  pattern 
depends  to  a  considerable  extent  upon  the  assumed  degree  of  the 
polynomial.  If  too  low  a  degree  is  selected  (e.g.,  linear, 
quadratic,  or  even  cubic),  the  resultant  deviation  values  will  be 
excessively  large  (thereby  requiring  a  less  sensitive  and  less 
useful  scale  factor  for  the  data  plotting).  On  the  other  hand, 
use  of  a  high  degree  polynomial  (e.g.,  9th  or  10th),  although 
capable  of  reducing  the  apparent  deviation  values,  probably 
detracts,  rather  than  assists,  in  the  data  interpretation 
process . 

The  accuracy  attainable  by  an  operational  ECXO  system  is  not 
necessarily  limited  by  the  apparent  deviation  relative  to  any 
polynomial  approximation.  With  a  suitably  designed  "look-up" 
compensation  table,  for  example,  the  individual  table  values  can 
more  closely  match  the  actual  frequency  behavior  of  the  resonator 
or  ECXO  oscillator.  The  ultimate  accuracy  will  then  be 
determined  by  any  lack  of  repeatability  in  oscillator 
performance,  together  with  the  quality  of  the  "look-up"  table  and 
associated  interpolation  procedure. 
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SUMMARY  AMD  CONCLUSIONS 
4.1  SUMMARY  AND  CONCLUSIONS 

The  program  has  achieved  a  major  objective  in  that  it  has 
produced  a  better  understanding  of  the  behavior  and 
inter-relationship  of  the  B-mode  and  C-mode  frequencies  of  SC-cut 
quartz  resonators.  The  investigation  has  demonstrated  that  high 
stability^  re-traceability/  and  accuracy  are  intrinsic  to  dual 
mode  operation  of  such  resonators.  The  experimental  results  show 
that  the  stability  and  retraceability  of  the  dual  B-mode  and 
C-mode  frequencies/  except  in  narrow  temperature  regions  where  a 
B-mode  activity  dip  is  encountered/  is  capable  of  providing  a 

g 

frequency  accuracy  of  several  parts  in  10  over  an  extended 
temperature  range  of  -55°C  to  +100°C.  This  potential  represents 
a  10  -  100  fold  improvement  over  conventional  TCXO  oscillators 
using  analog-voltage  varactors  for  frequency  compensation. 

The  B-mode  activity  dip  in  some  crystals  results  in  a 
moderately  small  frequency  perturbation  over  a  limited 
temperature  span  (e.g./  a  discontinuity  of  3  x  10  over  a  few 
degrees  Celsius).  In  other  crystals/  however,  the  interference 
may  be  so  severe  that  normal  oscillation  ceases.  The  question 
must  therefore  be  raised  whether  the  B-mode  activity  dip  problem 
can  be  eliminated  (or  circumvented)  by  modification  of  the 
crystal  design.  Is  it  feasible/  for  example,  to  design  an  SC-cut 
crystal  so  that  any  activity  dip,  if  present  at  all,  will  fall 
outside  the  normal  operational  temperature  range?  (Several  of 


the  resonators  used  in  this  study  had  a  single  activity  dip  ir. 
the  neighborhood  of  +65°C.  Could  this  interference  teiriperatur e 
region  be  increased  to  +85°C,  or  even  higher,  without  introducing 
another  dip  at  some  lower  temperature?)  Or,  is  it  feasible  to 
dimension  or  contour  the  crystal  to  suppress  interfering  modal 
signals  to  a  sufficiently  low  level  so  that  the  interference 
effect  will  not  be  too  great?  In  summary,  the  activity  dip 
problem  appears  to  be  the  most  severe  obstacle  to  the  use  of  dual 
B-mode  and  C-mode  frequencies  for  ECXO  compensation  purposes. 

On  the  positive  side,  the  B-C  differential  frequency  mode  of 
operation  appears  to  offer  several  significant  advantages  over 
alternative  techniques  of  measuring  resonator  temperature.  Use 
of  the  B-mode  frequency  provides  a  high  frequency  vs.  temperature 
sensitivity.  There  is  roughly  a  20  kHz  change  in  the  B-C 
frequency  over  the  -55°C  to  +95°C  range  with  a  5  HKz  SC-cut 
resonator  (corresponding  to  a  slope  in  excess  of  26  ppm/°C 
normalized  to  the  C-mode  frequency) .  This  high  sensitivity 
relaxes  the  resolution  and  precision  requirements  of  the  digital 
counter  required  within  an  operational  ECXO. 

The  excellent  performance  observed  at  the  lew  temperature 

extrem.e  was  not  anticipated  prior  to  the  investige  tion .  It  had 

been  thought  that  the  stability  and  accuracy  of  the  dual-mode 

compensation  technique  would  deteriorate  because  of  the  steep 

slope  of  the  C-mode  frequency  vs.  temperature  characteristic  in 

this  region  (greater  than  2  ppm^/C)  .  To  achieve  a  frequency 

—  8 

compensation  precision  of  ±1  x  10  with  this  severe  C-mode  slope 
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implies  a  thermometry  error  of  less  than  ±C.005°C  for  the  C-mcoe 
frequency. 

With  TCXO  oscillators,  the  crystal's  cut  angle  and  turnover 
temperature  m.ust  be  carefully  chosen  so  as  to  minimize  either  the 
maximum  slope  or  the  maximum  frequency  excursion  over  some 
selected  temperature  region.  From  the  results  of  this 
investigation  it  does  not  appear  necessary  to  put  the  same 
restriction  on  SC-cut  quartz  resonators  to  be  used  in  dual  mode 
ECXO  applications.  A  steeper  frequency  vs.  temperature  slope  may 
require  a  wider  range  or  higher  resolution  in  the  ECXO  look-up 
table  or  compensation  algorithm,  but  that  is  a  distinct  problem 
from  the  question  of  the  intrinsic  stability  of  the  quartz 
resonator.  Specification  of  the  lower  turnover  temperature  of  an 
SC-cut  crystal  intended  for  dual  mode  (B-  and  C-mode)  ECXO 
operation  will  be  determined,  not  by  an  intrinsic  instability  of 
the  quartz  material  at  either  temperature  extreme,  but  rather,  by 
the  minimum  spacing  of  the  stored  compensation  values  in  the 
look-up  table  together  with  the  quality  of  the  interpolation 
algorithm  used  for  intermediate  temperature  points. 

It  needs  to  be  recognized  that  the  factors  that  limit  t''e 
ultimate  performance  of  an  ECXO  oscillator,  including  the 
limitations  of  the  resonator  itself,  have  not  yet  been  fully 
established.  Design  criteria  that  have  been  developed  over  t:;.e 
past  decades  for  conventional  TCXO  oscillators  are  not  liKely  tc 
be  applicable  to  the  development  of  ECXO  systems. 

The  contract  study  has  also  demonstrated  that  thermal 
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transient  effects  cannot  be  neglected  if  high  clock  accuracy  is 
to  be  achieved.  Temperature  rates  as  low  as  ±5°C/hour  during 
temperature  cycling  produced  noticeable  frequency  offsets.  Sene 
of  this  deviation  may  be  due  to  "hysteresis";  however,  a  better 
method  of  defining  and  measurement  of  "hysteresis"  under  changin 
conditions  of  temperature  cycling  (with  variable  temperature 
limits  and  under  different  rates  of  change)  appears  desireable. 

The  test  results  to  date  show  a  significant  improvement  if 
some  degree  of  thermal  rate  compensation  is  used  to  augment  the 
fundamental  compensation  for  frequency-temperature  dependency. 
Further  investigation  will  be  required  to  determine  the  most 
appropriate  rate  compensation  algorithm  that  will  provide  "best" 
performance,  in  an  empirical  sense,  over  a  wide  range  of 
realistic  environmental  conditions. 

All  of  the  test  data  presented  in  this  report  are  shown  as 
frequency  deviations  from  either  a  5th  or  6th  degree  power 
polynomial.  This  procedure  allows  the  individual  raw  data  point 
to  be  graphically  displayed  to  a  high  resolution  level  on  a 
single  plot  covering  the  entire  temperature  range.  An 
operational  ECXO  unit,  however,  is  likely  to  employ  a  detailed 
look-up  table  in  which  the  stored  compensation  values  follow  the 
actual,  measured  data  points  rather  than  the  assured  polynomial 
approximation.  The  scatter  and  retraceability  of  the  data, 
rather  than  the  deviation  from  an  arbitrary  approxim.aticn ,  ther. 
provides  a  good  measure  of  frequency  accuracy. 

Inspection  of  all  data  for  the  various  runs  indicates  that 


an  average  frequency  accuracy  of  several  parts  in  10“  (for  the 
SC-cut  quartz  resonator  elements)  is  attainable  with  dual  mode 
frequency  com.pensation .  From  this  it  can  be  concluded  that  an 
ECXO  clock  error  of  several  millisecond/day ,  or  less,  is  possible 
if  the  resonator  performance  is  not  seriously  degraded  by  the 
oscillator  circuitry  or  external  compensation  processing.  There 
remains,  finally,  the  problem  of  aging  and  frequency 
recalibration.  All  operational  ECXO  clock  devices  must  include  a 
suitable  mechanism  for  both  periodic  time  synchronization  and 
periodic  frequency  calibration  (syntonization) .  Any  realistic 
estimate  of  clock  error  must  take  into  account  the  probable  (or 
worst  case)  error  in  syntonization  and  the  time  interval  between 
time  updating  under  anticipated  military  usage  conditions. 
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APPENDIX  A 


MEMORANDUM 

FROM;  Bob  Nelson 

TO;  Donnell,  Baltser,  Stone,  Quinones 

SUBJECT;  A  Hardware-type  ECXO  Digital  Compensation  Circuit 

Figure  A-1  shows  a  design  for  a  hardware-type  ECXO  compensation 
circuit. 

The  B-mode  (5.5  MHz)  and  C-mode  (5  MHz)  crystal  outputs  are 
first  prescaled  to  43  and  39  kHz  respectively.  The  39  kHz  is 
divided  by  the  lower-left  programmable  divider  to  1  Hz.  This 
divider  is  preset  with  a  value  between  26464  and  26480,  for  a 
division  ratio  between  39072  and  39056,  or  5001216  and  4999168 
overall,  respectively.  This  range  of  division  ratios  correspond 
to  the  range  of  C-mode  frequency  which  the  design  can  tolerate. 

The  exact  division  ratio  is  determined  by  the  value,  0  to 
16,  to  which  the  5  lowest  divider  bits  are  preset.  This  value  i 
the  5  upper  bits  of  the  21-bit  sum  of  a  20-bit  stored  correction 
with  a  16-bit  remainder  from  the  previous  such  sumiration.  The 
use  of  the  saved  remainder  to  bias  subsequent  corrections  n.eans 
that  the  16  lower  bits  of  the  correction  have  the  same  effect  on 
the  average  as  if  they  actually  preset  the  divider  with  a  4-fc:it 
integer  and  a  16-bit  fraction. 

The  relationship  between  the  stored  correction,  K,  and  the 
C-mode  frequency,  f,  may  be  seen  by  realizing  that  the  fixed 
preset,  plus  the  correction,  plus  the  counts  caused  by  the  39  kH 
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clock  must  total  to  the  full  range  of  the  divider,  2^®  Hz  or 
3.90625  X  10“^°  X  f  (at  f  «  5  KHz). 

The  sununation  is  done  in  serial  by  4-bit  words.  The 
impending  overflow  of  the  divider  is  detected  and  used  to  enable 
shifting  by  fours  of  the  16-bit  remainder  register  (lower-right) . 
As  each  hexadecimal  digit  is  shifted  off  the  right  end  of  the 
register,  a  correction  digit  from  the  27C128  memory  is  added  to 
it  and  the  sum  stored  in  the  left  end  of  the  register.  A  single 
54C74  flip-flop  saves  the  carry.  After  4  shifts,  the  register  is 
back  to  normal  and  has  been  incremented  by  the  16  lower 
correction  bits.  The  last  4  upper  correction  bits  are  added  to 
the  carry  in  the  last  period  before  divider  overflow  and  stored 
in  the  5  lower  divider  bits  along  with  the  constant  26464 
(6760^g)  in  the  11  upper  bits  when  the  divider  is  preset  at 
overflow. 

Only  part  of  the  16384  x  8  correction  memory  is  actually 
used.  The  corrections  are  stored  as  5  x  256  x  4,  but  another  256 
X  4  area  must  contain  zero  to  insure  that  the  carry  flip-flop  is 
cleared  before  the  first  summation.  The  memory  map  is 

ADDRESS  DATA 

oooooxxxxxxxxx  xxxxxxxx 

ooooioxxxxxxxx  xxxxoooo 

OOOOllAAAAAAAA  XXXXKKXK  (LSW) 

OOOIOOAAAAAAAA  XXXXRKKK 

OOOIOIAAAAAAAA  XXXXRKKK 

OOOllOAAAAAAAA  XXXXKKXK 

OOOlllAAAAAAAA  XXXXRKKK  (MSW) 

All  Other  Locations  XXXXXXXX 

A  ■  correction  address  K  ■  correction 


If  desired,  the  3  upper  address  bits  could  be  used  to  select 
among  8  sets  of  corrections  to  allow  a  small  ageing  adjustment, 
or  to  produce  intentional  frequency  offsets. 

The  corrections  are  selected  from  second  to  second  by 
counting  the  43  kHz  prescaled  B-mode  signal  in  the  8-bit  counter 
in  the  upper-right.  This  counter,  of  course,  overflows  many 
times  during  the  second,  but  its  final  values  cover  a  range  of 
about 

(5.5  MHz  X  2“^  X  1  sec)  x  (26  x  10"^/C  x  1400  s  156.  This 
assumes  a  26  ppm/C  frequency  vs  temperature  slope  for  B-mode  and 
a  140°  temperature  span. 

As  the  divider  approaches  overflow,  the  output  flip-flop  in 
the  upper-left  is  set  synchronous  with  the  39-kHz  clock.  Its 
output  is  then  synchronized  with  the  43-kEz  clock  and  the  rising 
edge  detected.  This  sychronously  resets  the  counter  and  loads 
its  final  count  into  the  8-bit  address  register.  Simultaneously, 
the  memory  is  enabled  through  the  end  of  the  summation  process. 
The  count  copied  into  the  address  register  is  always  one  less 
than  the  correct  value  (module  2^^},  because  one  count  is  lost  to 
resetting  the  counter. 
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ZZZ  PRINT  "t  tOPERATOR  MUST  INPUT  SYNTHESISER  FREQ  'NOMF'.:;)'  AS  STORED  IN  *" 

Sf.4  ®PINT  "t  RAW  DATA  FILE  ’  HEADER!’  (SEE  PROGRAM  •' DATA/PLOTl  ’  ;  .  ALSO.  *  ' 

PR  IN^  "t  OPERATOR  MUST  INPUT  INFO  ON  CRYSTAL.  SELECTED  DATA  PAIRS.  ETC  «  " 

0_-0  PRINT  i[*tt*«*t»»»«»«*»»»***»»«»*»*»«********»»*»*»*^»****»*P*******'«»T.  ■ 

_Tv  PRINT 

TOV  DI3P  "INPUT  OISK.Ws  ?  .DATA  FILE  NAME:  ?"  9  INPUT  DSK*.DATFIL» 

-M-i  'sc-lMT  "DATA  STORED  ON  DISK:  “;DSK*;TAB  (40)  ;  DATFIL#:  TAB  (60):  DATE  :TIME 
-'sv  DIM  A(  1,')  .R(  .'.GJ  .  r<3)  ,B(3100,6)  .CaMMeNTI*C20OJ.C0MMENT0*CS00I.COMMENTr*C;'>;  : 

ZZ'j  REM  *  f  Above  DIM  A.R,S<T  permits  up  to  6th  degree  pol  vnom:  al :  D IM  B  permits 
.'I'.’’  lines  o-^  raw  data  *  6  variables  (inc  TIME) 

74;  CISP  ■'  INPUT:  CRYSTAL  TYPE:. SERIAL  NO.  FREQUENCIES.  ETC 13  CHAR  MAX  COMMA)'" 
INPUT  T'^PE*.  serials.  FREQS, OVRTONEi  .  OTES,  years.  ETC* 
print  ■■  crystal  TYPE: 

'  ■•  •  PR  1  NT  TYRE* .  3ER I  AL  * ,  FREQ* ,  OVRTONE*.  DTE* .  YEAR* .  ETC* 

:0"  '^OF-  J»1  ro  4 

DiSP  "INPUT  SYNTHESIZER  FREQUENCY.  NOMF  (  "  :  J ;  "  )  :  "  '?  INPUT  NCMF  '  J  ■ 

-t-lMT  "  S/NTHESIZER  frequency.  NOMFC’:  J;  ••)  :  "iNCMFIJ: 

'iE,>r  J 

:_V  L'lrP  "ANY  ADDITIONAL  COMMENTS  (UP  TO  COO  CHARACTERS — NO  COMMAS'.'" 

471  I, -PUT  COMMENT!  » 

i-  PiNT  "NOTE:  ";  COMMENT!* 

15"  PRINT 

J-  CISP  "DEGREE  OF  EQUATION''" 

4-"  PRINT  "DEGREE  CP  EQUATION": 

IS"  INPUT  D 
4.1,.  pp’INT  D 

5  ;'  CISP  "TIJTAL  NUMBER'  OF  RAW  DATA  POINTS  IN  RUN'""  '  Must  De  1  mown  zr  esTimalsc 
51U  PRINT  "TOTAL  NUMBER  OF  DATA  PAIRS*": 

-IC'/  input  total 

-:z:’  PRINT  TOTAL 

"■4.,)  C'ljp  "SEblNNING  DATA  PAIR  NO.’f-for  deriving  oolvnomiai  tset  -  i  ci  .=n  t  a  • 

55:'  iriFUT  .. 

5i"  PRirir  "BEGINNING  DATA  PAIR.L*":L 

=  o;,  ip  "Si'iDING  DATA  PAIR  NO. "''.for  deriving  BOlvnomial  coer  f :  c  i  in  t  s ' 


pc'iMT  "ENDING  Data  pair.  M*";m 
'  OISP  "STEP  3I7E"''/Gr  deriving  oolvnomiai  coe-f  t  i  ■:  i  ant  3 )  " 

_  1  ■'  I NP  U  T  5 

•;C  .iciriT  'STEP  •S::E.5*";S 
,r  .1=1  *!Nr  I  .  M-L.'  '5) 

.J,,  c.i;N:r  ’MUMBEP  IP  POIMT'j  '^SEC  '0  DETERMINE  POLYNOMIAL  =•  I  ^  r 
.5'  OI.ip  'ii.iMEEP  OF  POINTS  USED  '0  DETERMINE  POL'-TICM  I  AL  ■^:"='t.N 


B-L 


MQTEsMust  match  disk  ^ile  namel  ic****t**'sitii* 


;.o  :■  F=:r:  iT 

"NO  ;  frtt  ■,  7  ■  :  “  T -FREQ H  Z  > '• :  TAB  Qi  ;  “F-FREa HZ  )  "  ;  Zi  '  :  "  ■■I-: 

!a  :  "  V, HI;";  tab  ;  5~  '  :  "  7  i  I1£  "  :  TAB  ;  Sdi  ;  "TEMP  >  APR  ROX  ) '• 
r'  H';-  A  ■  i  I  *N  ?  r  D*  2: 

PCP  TO  Z«D+1 

A'J'=0 
’1-:  .•■(£>  T  J 
Z-.-'  FOR  K=i  TO  D+l 
^  T  I'  )  =o 
“ -i'j  NEXT  K 

'T'X  ASalONlt  1  TO  “ XTALDAT;;  '  '  MQTEsMust  match  disk  t^ile  namel  itttttttrrt* 

FOR  1*1  TO  TOTAL 
-X  FOR  J*l  TO  -a 
’iO  READ#  1  :  Bd.J) 

■  NEXT  J 

e -  r .  i  > »ncmf ( i  j  / 10.24*  (B( i,  i)+iO  '7)  i?  ad. 3) •nomft  1  >/ 10.24*  o d.  z)  -mo'  7 > 

31'/  B  (  1 . 2)  “NOMF  (2) /lO.  24*  (B  d  ,  2)  ♦lO''?)  !  Translates  to  F-mode  ■frequency 

32;'  S  d  .  3)  »NOMF  <  4)  /  10. 24*  <B  d  .  3) -klO^'?)  '  Translates  back  to  C-mode  -frea 

33'.'  3  d  .  4;  *NOnF  (3)  /  10. 24*  <B  (  1 . 4) +10-'‘7)  !  Translates  back  to  B-mofle  •freq 

.1 4-;  next  I 

330  “Qp  TO  M  STEP  S 

da'.'  <=  I  S  ( 1 .  1 )  ♦e  ( 1 , 3j  ) /2-S  d  ,  2) -XCONST  !  Operator  must  replace  XCONST  with  a 
numeric  chosen  to  reduce  magnitude  o-f  X  over  the  temperature  raange 
'3~0  v*6 d , 2) -TCONST  !  Operatore  must  replace  YCONST  with  a  bias  so  as  to  reduce 
■'  tc  near  0 

23'  ’  TAP^ROX*- (B  d  .  4) -B  d  ,5) -470000) /13a)+aS  !  Uses  old  B-C  mode  temp  aporo;; 

23':  PRINT  USING  900  )  I .  ^B  (1 .  1 )  ■►B  (1 , 3)  ) /2,  B  (1 . 2)  ,  X ,  Y.  B  <  1 . 6)  .  TAPPPOX 
['"AGE  40.2X.3D.4D.  IX, 30.40.  1X.S4Z.3D,  1  X.  S3Z.  40.  3X,  50.  dX  .  S3D.  20 
-  L'.'  REM  t********************* **«**««**«****««*«« t************************** 

-20  REM  «  FOLLOWING  PORTION  USED  TO  DERIVE  THE  POLYNOMIAL  APPROX IMATICN. P  r 
REM  **************************************************** ******»«**«tit*it* 

■"4'>  F'JP  J*2  '■Q  2*0+1 
A  .  J  '  »A  ( J  )  +X-''  ( J-1 ) 

'4,"  I'.E-'.r  J 

■  for  h*l  TO  0+1 

■S'!'  P  '.h  .  0+2)  "T  '.K  )  +Y«X-'  (K-l ) 

-ip-.'  r  '.K  )  -T  (,k  )  *ytx  "  '.K-l ) 

NE'.T  k 

r:';-;.'  r d+2)  »t  <0+2) +y''2 

..■2'-'  NEXT  I 

I.’:-;  FOR  j*i  ro  o+i 

:  '40  FOR  K«i  ro  0+1 

R  J  .  K)  -A^  J+K.-l  ) 

L-'O'.''  NEXT  H 

next  j 

I  FOR  J*1  TO  0+1 

1  !•  aj 

:  ’  '.''I'  IF  R  ' ,  J  )  . ..  0  THEN  1  1 50 

I  1  i'.i  )■.*•<  +1 

:i2'j  IF  .  0+1  THEN  1100 

.170  PRINT  "  NO  UNIQUE  SOLUTION" 

114'.'  GOTO  17*0 
1130  -R  1*1  TO  0+2 

I I  ok'  G*F; J .  I ) 

1  170  R'd,  I  *ft'l<  .  I) 

1  1 3i.'  R  ',  k  ,  [  '  *S 
.1-^0  NEXT  I 
1  _  .'<•>  2  *  1  -  R  J  .  J  : 

I- -10  FOP  1*1  ro  0+2 
•  =  ■■  J.  1 1 »:  »p  ; 


I'lExr  1 

,-OR  f-=l  TO  U*l 

rnsr’i  1300 
R f  . .  J ) 

:3~0  FOP  1=1  TO  D*2 
:  ZS’J  F.  ■:  K  .  I  I  =R  ■  K  .  I )  +Z  «R  J .  I ) 

IZ~0  ME XT  I 
trov  NEXT  K 
iTxo  NEXT  J 
.320  PRINT 

•  :T0  print  '•  CONSTANT  (AO)  «'* :  R  <  1 , 0+2) 

-.340  FOR  rO  D 

1 :5o  PRINT  J; “DEGREE  COEFFICIENT  (A":J:“)  ; R ( J-d  . D+2) 

1  ;aO  NEXT  J 

PRINT  "  'Y (IDEAU)=A0+<A1*X)+<A2*(X*X) )+ . (An*(Xn))  “ 

1  •30  PRINT 

I  39i_i  o»(j 

:  400  FOR  J=2  TO  D-t-l 

1 4i0  F»P+R ( J.D+2) « (T ( J) -A ( J) *T ( 1 ) /N) 

142'.-  NEXT  J 

1430  Q»r(D+2)-T<l)  '2/N 
1440  Z-G-P 
1430  I»N-D-1 
144jO  PRINT 
1470  J=P/Q 

■430  PRINT  “COEFFICIENT  OF  DETERMINATION  (R'"2)=";J 
1  ;P0  PRINT  "COEFFICIENT  OF  CORRELATION  -••;SQR  (J) 

IS'XO  PRINT  “STANDARD  ERROR  OF  ESTIMATE  (HZ)  ■••;SCR  (Z/I) 

.SI  ?  PRINT  “STANDARD  DEVIATION  C'F/F)  ■‘•;SQR  ( Z/N)  XNOMF  ' 2)  /  10-  6 

IS20  RE.M  *<*****END  OF  POLYNOMIAL  ROUTINE  **»««t*»»**1t*»**«***t«*.t»***it*f * 
1330  PRINT  (J  PRINT 
tS4<;i  SEEP  30.2000 
1330  ALPHA 

'So*?  DISP  "INPUT  NEW  3E3INNING.  ENDING  DmTA  PAIRS.  AND  STEP  SIZE  OR  •RESET 
MC  PUN)" 

13“'?  INPUT  L.M.3 

1330  PRINT  “  BEGINNING. ENDING  DATA  PAIRS,  AND  STEP  SIZE:  “:L:M:S 
iS'O  PRINT  I*  PRINT 
ic*.*'?  BEEP  1'.’*.  1(1*0 

•ai;*  DISP  "INPUT  DATE  OF  RAW  DATA  RUN  (IN  FORMAT  ’DAY  MONTH  YEAR')" 
liC?  INPUT  PDATE* 

la:'?  orsp  "ANY  ADDITIONAL  COMMENTS  (UP  TO  200  CHARACTERS — NO  COMMAS')" 
lo4(,’  input  CCMMEMT2S 

LiS'.'  PRINT  "  NOTE:  ";CaMMeNT2»  U  PRINT 

•.=a;i  PRINT  TAB  (2)  ;  "NO.  "  :  TAB  (7)  :  “  T-FREQ  (HZ )  " ;  TAB  (20)  :  “F-FREO  '  HZ .  ' :  T^e  3 
«;■  "rrAB  (43'.  ;  "C  (  IDEAL)  ";  TAB  (33);"  DEV  (  -F/F)  "  ;  TAB  ■  66)  ;  “TEMP"  :  TAB  (  73  ■  :  ’ 
;r7-;.  PRINT  " - 


.  od'.'  3='?  H=(;'  S  PRNT*0 

la'*'."'  FOR  I=L  TO  M  STEP  3 

1-.,..^  Na;i-INT  ((M-D  -'S)  i>  PRNT*PRNT*t 

;  7  (■,'  (=  '  B  '  1 .  1 -yD  '  1 . 3)  )  .  2-B  ( 1 . 2) -XCONST  '  DiP  +  erancs  fraa.wit.M  jT4.=  at..n  mc 
1'-'-’  '( =B  (  1 . 2 )  -  YCONST  I  'CONST  to  6a  raolacad  with  numaric 

1  Ti.'  TapPPOX'*- (  '. B  ( 1 . 4)  -B  ( 1 , 3) -47(;)'IiOO)  / 1 33) -ySC  '  TamDeratura  aopro;;  casec  4i 
Jift-eranca  Prequenc, 

;  ■  1...  p3iP:  <  1 , 0+2) 

1  3*.  ruP  >J  =  I  TO  D 

~  -'^P  *P  .  J  + 1  .  D+2 )  *  X  '  J  '  Simple  tamoaraturs  appro:,  (wi  CTiOut  thernal  rata 

NEXT  J 


2  '!•  G=G+ABS  ;  Y  -P ' 


;'ii>p  usiNii  ;  i .  o  .  i -►e  3)  > /2,S’.  i.  •  mci^p  _  ■ .  i 

iz'.'.'  tF  ’^RNr^l'j  THEN  1310  EL3E  1340)  '  Prints  avsrv  lOit."  iii-ia  in  same  1  a staa  S' 
:zl.  PFINT  CEING  1320'  ;  I .  '  B  <  1 ,  1 )  >8  <  1 . 3)  )  /  2,  S  ( 1 . 2  J  ,  .  P.  Y-=' 'NCITF  ■:  2 )  .  1  •  =.  TwpPF.C 

a-''-  InmGE  4D.  IX.  3D.  30.  1  X.  8D.  3D.  l.X.  34D.  3D,  IX.  34D,  3D.  1  <.  3 1  Z  ,  2D9 . 2X  .  332 . 20 .  O'  .  iO 

.3  30'  ■=F'Mr=0i 
134"  -JEX  !'  I 

'25y  'J=')-'.'N)  '.S.NOMF'.  2)  '10)-'6  i?  F*Q,'N/NOhF  i2) /' 10''6  iB  PRINT 
.ja'j  PRINT  TAB  '.SJ;'"***  standard  DEVIATION  (''F/F)  »  ";U 

a’‘.)  PPXMT  TAB  AVE  ABSOLUTE  DEVIATION  ('F'F)*  "  s  F  R  PRINT 

la-ao'  aesuB  2i  =u 
;a°'.  jGEUB  1930 

BEEP  30.2300  9  GOTO  1330 
PRINT  '■ 


L=2:!  END 


19-.;)  ..;0LD»O 

1940  FOR  I^L  TO  M  STEP  S 
>  Q20  •<=  'B  ( I.  1)  1 ,3)  J  /2-B  (1,21  -172000 

L  R  ia-j  y =B  (  1 . 2 )  - 1 000 1  OOiO 

19-j  TAPPROX*- ( (B < I . 4 J-B ( I , 3) -470000) / 133) >83  !  Tamearatur*  appre:;  basaP  cn  B-C 
39  sacand  SC-cut  «taX 
1930  P»R(1.D>2) 

1- 590  FOR  j-l  TO  D 

ZOOO  P=P-*-R  ( J>1 ,  D>2)  «X-J 
2010  NEXT  .J 

2'.' 20  '-■»  <  r-P>  « lOO/NOMF  (2)  I  Sata  vartical  scale  to  3»10''-S 
y.'ZO  PEM  Oelatad  ■'  PRINT  IjX.V.P  ’  on  2  Dec  93 

nave  <.v 

_-:-S0  IF  XCLD'-  X  THEN  2060  ELSE  2070 

2':- 50  L.aBEL  ■i  BEEP  10.  lOO'  '?  GOTO  2080  !  Labels  decraasinq  tamp  soints 

Z'l'TO  L.43EL  "o"  iJ  BEEP  13.100  I  Labels  incraasinq  tamp  points 

2- .)£0  'CLD=X 
2-.-90  mE.XT  I 

_  1 '.)'!'  BEEP  40 .  1 000  i*  BEEP  30 ,  1 000 
21  WAIT  20000 
2120  DUMP  GRAPHICS 
2130  ALPHA 

214<.)  OISP  "ADD  COMMENTS  (REPEAT  TYPE  OF  CRYSTAL,  DATE  RAW  DATA.ETC(NC  COMMAS;  ' 

2 150  INPUT  COMMENT'S 

2130  PRINT  I?  PRINT  C0MMeNT3S 

2170  RETURN 

213'.'  REM  **«*******. «*********»**»»*»****»**»*.«»)([»*»«»«********. -4* 

21^0  PEM  *»rt  FOLLOWING  PLOTS  THE  FRAME  AND  AXES  *  ■« » 

I'.K  '  REM  t*«*»*««»*»»***»****»*»»'»*»*********«****»**»*4***»**,«'l*ttT*-r««4«-» 

221':'  CLEAR 

-220  P').JTTER  IS  1 

223-.I  2-=3  )t4«'«4t«C*3  qivss  -full  scale  vartical  Y-AXIS»-^/-3-*l'.')  -5 

—2  4'.'  ijPAPHALL 

223'.  uCCATE  O'.'i .  O'JO .  1 3 , 30 

22i5C  -prtME 

227':'  DEG 

223'.-  PEM  X*  <MAX  and  .'.MIN  must  ba  modi-fiad  to  match  observad  .'.MhX  and  'MIN 
22-?'.  PEM  *t  «iso.  all  other  horicontal  scala  values  mav  iieed  to  ba  Tiopined  tt 
■'i.e  dcrract  label  positioninq 

-3  .)'.)  »Max  =  1'.'00  I?  XMIN=-1200  '  Ad  lust  to  observed  XMAX  and  xniM  »-'dr  test  .  tai 
231  '  CSI2E  i 

2320  SC.'iLE  'MAX  .  XMIN.  -C.  C 
233.'  'AXIS  'MIN.  .  5. -C.  C 


-r'S'.  .  '  1 . 45 hC) 

Cil-H  4 

Li'SEL  ■'■JFESET  r-F  FREQUENCY  vH2'- 
NO ve  X  N(i  X  - :  4i;io .  -(.3*0 
0  3 1  »  »  4 

O.NtiEl.  ''RAW  OAFAj  " :  RDATE’i 
NO'v'E  <MMX-600,  1 . 4«C 
C3IZE  4 

LAEEL  "APPROX  rENPERATURE  tC) " 

NO'-'E  X  N A  X  -►  1  O' j .  -  ( .  93  «  C 
CEG 

LOIR  90 

UAbEL  '  F/F  FREQ.DEVIATICN  (10  )  " 

MOVE  XMAX  +  113, -PSliC 
LABEL  ’-3" 

MOVE  XMAX+1.0 
LINE  TYPE  3 
DRAW  XMIN.O 
T-=-33 

FOR  I-*0  TO  31  STEP  1  '  *  Qiva*  temper aturs  span  from 

X=*'.MHX-l';i-I  «  (  XMAX-XMIN)  /  31 

MOVE  X . 1 . OS*C 

LABEL  T>S#I 

MOVE  X. ,93*0 

(.ABEL  "-" 

NEXT  I 
RETURN 


APPENDIX  C 


3^S44  .  “T3 

iriT  »»:«.»«****?■»«*»****» 

"■«  PRaSRAn  •  DATA/PLOT:  '  revised  ’-10-S4 

:  -'C- '-IT  "It**  **»««****»»»  :*»***»»■»»*««****»****«*«(! 

:r:'  =RiMT 

!.  (0  RRIMTER  is  '.’03 

15'.’  i^RINT  "******:« ***************************************** *****************<*' 
liV  PRI.MT  THIS  PROGRAM  ’  DATA/PLOT::  '  CONTROLS  4  CHANNELS  OF  MODEL  DIO”  *“ 
PRINT  "*  REFLECTOMETER,  STORES  a  RAW  COUNTER  READINGS  UN  DISK,  PRINTS  GUT  t  ■ 
13V  CRINT  B?<C  MODE  FREQ  ( S.  4a?<3.  OOMHZ  IN  THIS  RUN)  FOR  THE  TEST  RESONATOR  « " 
1-'.’  PRINT  "*  (CHANNELS  0  S<  1)  AND  THE  B  ?<  C-MODE  FREQUENCIES  'OF  A  REFERENCE  *  " 
2'y.'  PRINT  •'*  CRYSTAL  (COLORADO  CRYSTAL  CO  ;S/N  1)  USED  TO  MEASURE  APPROXIMATE  T  " 
GI-.'  PRiNf  "t  OVEN  TEMPERATURE.  THE  T-MODE  FREQUENCY  OF  THE  TEST  CRYSTAL  IS  *  " 
_D  .  P'^  INT  "t  assigned  to  CHANNEL  Ot  F-MQDE  FREQ  TO  CHANNEL  1.  THE  B-MODE  *  " 

J3  ’  -'PINT  "«  OF  THE  REFERENCE  CRYSTAL  IS  ASSIGNED  TO  CHANNEL  D;  C-MODE  ''0  »  ' 

PRINT  '•*  CHANNEL  3.  IF  THE  PROPER  POLYNOMIAL  COEFFICIENTS  ARE  KMCWN  FOR  *" 
PRINT  "»  THE  TEST  CRYSTAL  AND  KEYED  IN  BY  THE  OPERATOR.  THE  FREQUENCY  *" 
Ci'.'  PRINT  •■*  DEVIATION  OF  F-MODE  FREQ  RELATIVE  TO  THE  PREDICTED  POLYNOMIAL  «" 
PRINT  "t  VALUE  'WILL  BE  PRINTED  OUT  AND  ALSO  PLOTTED'.ON  THE  HP-ST  CRT)  *" 

_3C  PRINT  “**»********»*»»*»****»*»**»«*«*»»»*»»»»*.*»»*****«**»*»»********»***" 

CP'.’  PRINT 

30'J  PRINT  " tit*t*t*t*»*tttt*tt***X**X*t**tt*$***tMt**ttX*t**ttt**ttttttt*t*tt*" 

no  PRINT  "t  NOTE:  VARIOUS  CONSTANTS  IN  THIS  PROGRAM  MUST  BE  MODIFIED  TO  ^IT  * " 
TDj  PRINT  '•»  THE  SPECIFIC  CRYSTALS  TO  BE  USED.  CONSTANTS  USED  HERE  APPLY  TG  t" 
PPINT  "*  COLORADO  CRYSTAL  S/N  2  IN  CHAN  0/1  COLORADO  CRYSTAL  5/ N  1  IN  *" 
740  PRINT  "*  CHAN  2/3.  ALSO.  THE  OPERATOR  SHOULD  CAREFULLY  FOLLOW  ' INPUT'  «" 
750  PRINT  "t  INSTRUCTION  ON  CRYSTAL  I . D. . SYNTHESI ZER  FREQUENCIES (MODEL  CIO”).*" 
74'.'  PPINT  PLANNED  TEMPE, -.ATURE  PROFILE.  ETC.  THIS  INFO  TO  BE  STORED  ON  FILE  « " 
7”  ;  -PINT  NAMED  ’  HEADER:: ALL  COUNTER  FREQ  READINGS  STORED  ON  '  XTAL/DATA;;  *  " 

=  '■'  PPINT  •'«  FILE.  i-l';i-34  t" 

PRINT  ••«**«*****»»««**»»*******************»*****»**«#********«, ******* 7***" 
."0  PRINT 
4  i  ■.!  PR  I N  r 

:D'.'  DIM  rPRCFILESClAOI,  IDENT1  (4)  C  1283  .  MODE*  <4 )  Ci43  .  COMMENT!  SC  160 3  .  R  :0.  lO) 

-''0  MASS  STORAGE  IS  "iDTOij"  !  Uses  floppy;  Winchaster  repuires  diffarsnt  addrsss 

45'.'  CREATE  "HEADERI  :  0700“ .  1 ,  1023  !  *CMang«  to  ’HEADER'/’  or  purqa  fila  if  HEADER! 

.slraadv  usaO  as  data  fila  nama 

46'!  ASSIGM4*  1  TO  “HEADERI"  '  *Must  match  nama  in  abova  line 
'  ■  >  CREATE  “  XTALDATA;  DT'l'O"  ,  7,  aOOO  I  Chack  CAT  to  ansura  no  duolication 
430  ASSIGN#  2  TO  " XTALDATA; 0700"  '  Ooan  a  floppv  disc  fila 

aP"  oisp  "Rftw  DATA  TO  BE  STORED  ON  D ISK#::  AND  FILE  NAME"';"  i  INP'JT  FLOPS 

:i'  INPUT  PHILE* 

C;".'  PRINT  I?  PRINT  “  DISK#";  FLOP*;  “FILE;  ";PHILE*;TAB  (65):DATE  ;TIME 
5:'.)  DCSP  “INPUT  PRESENT  DATE  (IN  FORM  '21  DEC  34’)"  '?  INPUT  DTE* 

5D';'  DI3P  "INPUT  APPROX  TIME  THAT  RUN  WILL  START  (  INCLUDING  ’AM’  OF:  PM')  '  INPUT 
"'ME* 

•77'.’  DISP  “PLANNED  MAX  TEMP  (Number  on  1 '/-—(daqraasC)  “  's*  INPUT  TMAX 

54'.'  Di3P  “PLANNED  MIN  TEMP  (Number  onlv — daqraasC)  “  'S'  INPUT  TMIN 

550  -PINT  "TMAXa  “;rMAX;“  aC“:TAB  <20);“TMIM»  '•;TMIN;"  oC" 

54';'  DISP  “PLANNED  TEMERATURE  PR0FILE-(2  tvped  lines  ma;:  '  “  '?  INP'JT  TFRCFlLEs 

5”-:  i:'PINT  “FLOPPY  D I SK# ";  FLOP'S ;  TAB  (43);  “PUN  STARTED :“;  DTE “  ":”''  MES 
53'  P^-'iriT  “PLANNED  TEMPERATURE  PROFILE;  “ ;  TFROFILE* 

5'’'.'  0 ISR  "ADDITIONAL  COMMENTS'’  (Ub  to  C  tvoad  lines;  if  nora.kavin  '’I’lno  ddiiiiTias 
■  ■'  ;■  INPUT  COMMENT ;S 

PPINT  “NOTE:  “;C0MMENT:*  ij'  DISP  '?  DIS.^ 

4i.'  PPIN'''#  1  ;  FLOPS.  PHILES.  DTE*.  TFROFILE*.  T't'MES.  CCMMENTl  S 

rC"  pi^INT  “CHArJ“:TAB  '.  7  '  ; '"CR'/STAL  MFRI-IDENT"  ;  TAB  ',  SO  '  :  '  MODE  "  :  TaE  '  NCM  .  =•  n  7:- 

rr.  4Q 


C-1 


I 


>"=’  r  ’  ri^sumds  ar  I  /  4  r5+ 1 -actcmetar  channels  in  use.  .-<1  ;n  ^.“linnsi 

I'.jnn^i  ;  UH3C!  Cwics  in  each  c/ule  ta  crov.de  i  -aw  rjata  channels 
■  ;:3F  ■IHfii'lMELW :  : :  "  :  INPUT  CPVSTilL  MFRi<  ICENT  'But  no  Commas''"  i-  INPUT  luEN 

'  ai  ,  ' 

■I-.'  I  ISP  "CH>4NMELtf :  I  ;  :  INPUT  OPERATING  NODE  ?;  OVERTONE/ FUND  (13  onar  /na.:;":!'  : 

;:-r  'IQCES  ,  I. 

-  =  '/  JlSP  "CHANNELS":  I :  :  INPUT  ASSOCIATED  SYNTHESIZER  FREQ  (In  MHz  units - lo  alo 

hi  inaractsrs  :  )  "  :  I"  INPUT  NQMF  (  I ) 

u  •  ■  '-RINT  TAB  I;  TAB  ( 7)  ;  IDENT»  ( I )  :  TAB  (45)  :  MODE*  ( I )  :  TAB  ~0 )  ;  NQMF  i  I  j 

u3';  '-Pints  l  ;  I dents  '  I  >.  mode*  <  I ),  NONF  ( I 

jj-.  ^„£yT  I 

■  -OP  lac  I  TO  3 

CIoP  -tittzs*  BE  SURE  TO  KEY  IN  SOME  ’NUMERIC’  FOR  THE  FOLLOWING  zrttftf 
•Uv  uISP  "INCIDENT  CRYSTAL  POWER  (MICROWATTS)  :  CHAM#":  I  INPUT  PWR I  ) 

■  :■/  PRINT  TAB  1(5)  :  "INCIDENT  CRYSTAL  POWER:  CHANNEL#"  :  I ;  PWR  (  I  )  ;  ‘MI  CRCWa  TT  " 
•’4/  print#  1  :  PWR(I) 

■5'.)  NEXT  I 
“be.  PRINT 

print  ASSUMED  POLYNOMIAL  COEFFICIENTS  »**«*» «***«*» it **** 1 1 ** « t t * » 

(tttt?  " 

”30  axioms  far  a  6th  ord#r  paXynomial 

•90  Disp  Hjp  pQj_Y  COEFFICIENTS  ARE  NOT  KNOWN.  KEY  IN  ’0'  FOR  EACH  QFTHE  FCLLOWIN 
6;  " 

aov  DISP  "  INPUT  THE  DESIRED  POLYNOMIAL  COEFFICIENTS  .  THE  PROGRAM  '■ 

3 15  DISP  "  ALLOWS  SIXTH-ORDER  COEFFICIENTS  TO  BE  KEYED  IN. 

CC'5  DISP  "  FOR  A  FIFTH  ORDER  POLYNOMIAL  .  KEY  IN  A(6)a0."  i?  DISP  '?  DISP 
33'.'  DISP 


‘J40 

OISP 

"CONSTANT 

A  ( 0 ) : " :  1? 

INPUT 

R( 1.0+2)® 

PRINT 

A(0)'a":R'l,D*2' 

’iS'l* 

OISP 

"COEFF. 

A 1 ) :  " :  ij 

INPUT 

R<2.0+2)® 

PRINT 

A(1)»":R(2.D+C) 

r™^ 

0 1 SP 

"CQEFF. 

A  (2.1  :  " ; 

INPUT 

R(3.D+2)® 

PRINT 

A (2) ; R (3. D+2) 

DISP 

"COEFF. 

A (3)  :  ";iJ 

INPUT 

R(4.D+2)® 

PRINT 

A(3)=":R(4.0-C) 

.59’  ’ 

DISP 

"COEFF. 

A  ( 4 )  :  " ;  1? 

INPUT 

R<5,D+2)® 

PRINT 

A(4)»";R':5.D+C; 

**•  •' 

;SO.  • 

DISP 

"CQEFF. 

A\5)  :  ":i? 

INPUT 

R (6,0+2)® 

PRINT 

A(3)»":R(6.D*2:' 

•  4  r'- '  } 

OISP 

"COEFF. 

A  6 )  :  "  :  1? 

INPUT 

R(7.D+2)® 

PRINT 

A(6) =" : R .7.3*2) 

^  1 '  ■  uT'  !<’E 
■I'/  FOR  I 
^■*'5  PRINT 
•4u  NCX ' 
'5"  PRINT 


YU  7  GO 
a  I  TO  7 
#  I  :  R 
I 

"tZZZZ 


GOTO  S4<5  li  WAIT  10000  !  USED  TO  ’ESCAPE’  ANY  TYPING  ERROR 

/ 

Rc.o+i:, 

ZZ  F=  A(0>+  A(l)»X  +A(3)«X'2  + . +A(ni*X'n 

TO  1 


h’5  ASSIGN#  1  TO  X 
PRINT  I?  PRINT 

Sv  PRINT  "MO.":TAB  (7)  :  "T-MCDE  (HZ)  "  ;  TAB  (2C)  ;  "''T"  ;  TAB  (03)  :  "F-MODE HZ  ■  "  :  ~AB 
:'F,L£F";TAB  ( 4= ;  ;  "TAPPROX  "  ;  TAB  ( 60 )  ;  "TIME"  :  TAB  (70)  :  "-'F/F" 

’0  0:3P  "THE  FOLLOWING  PORTION  INPUTS  VARIOUS  FREQUENCIES  FROM  ’’HE  HF533uA ' 

'/■;'/  OiSP  " COUNTER  AND  STORES  DATA  ON  THE  FLOPPY  DISC.  A  DATA  FILE  CAuLEZ  " 
.1'.’  OIuP  "  *(TALDhTD*  is  CREATED,  check  CAT  -for  duel  i  cati  on  '  ' 

01'.<  GOSUB  1670 

'-'.I  LIrlE=>.i  'Y  iK:iP=0  R  X0LD*'5  i?  C0LD*0 

'4.;.  Fe,i>i  z  zzztzzzzzzztzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzztztzzzztzsz  TIHTtHTtllc  ft 
jS'.'  =:EM  ttz  FOLLOWING  SUBPROGRAM  CONTROLS  THE  CIO"’  PEFLECTOMETER.  OCUNTEP.  s  r -- 
REM  «.t«1i**'«******T*«*««1[*»«*»«t»1i***,tt**t*«**ii****it*«it»f»»?itt«((«*»t*«t» 

UN  riMEF:#  1 ,  S'OuOO  GOTO  I'jPu  '  Gives  50  second  cvei  a  tiine;  na-/  oe  mcci-tac 
"3'.'  GOTO  1030  '  Idle  loop  for  timer 

OUTPUT  XOfii  :  ■93"  I?  OUTPUT  7'5o  :  "HH"  ’  Mul  ti  p  1  a:t  an  selects  CHANNELS  313 
1  ."j  OUTPUT  ■’■jh  :"FN1"  '  Switches  to  COUNTER  A  channel 
1.'  I4A !  r  '  Reduced  tram  earlier  3<500  value«t 

liNTER  :  CC  '  Counter  readme:  Por  Channel  C  to  computar 

:  r ulJTPUt  ■’'5'  ;  ■'•'Ml'’''  '  Switches  to  COUNTER  B  channel 

;  LI'.'TER  T;.!-’  .  '  Ccunter  readme  tor  Channel  3  tc  cemputer 


I 

1  O’." ’='.17'  7.;-7  .-"'-Ml"  ■  Switch  back  tc  COUNTER  .i  channel 

ENTEF  ~'.'i  ”  :  Ci-iC  '  Cnd  raaaina  o+  Channel  2  sent  to  conouter 
ll.?-,'  OL.TFL'T  ;  '^3"  !  MultlBleiisr  selects  CHANNEL  O  M 

; :  - 'jurPLiT  7'j7  ;"Fmi"  '  awitches  to  COUNTER  A  channel 
.  C  :v;-iiT  SCO  '  Reduced  tram  earlier  SuOO  value 
I  EriTEF;  ’0~  ;  co  ‘  Counter  readinq  tar  Channel  0  to  camoutar 

ICC"  OUTPUT  7'j7  ;''FN17'‘  I  Switches  to  COUNTER  B  channel 
'  C3'  ENTER  TO’  ;  Cl  I  Counter  readinq  -for  Channel  1  to  comouter 
iCP'O  OUTPUT  707  ;"FN1''  I  Second  switch  to  COUNTER;  A  channel 
;CS'.‘  ENTER  707  :  CaO  '  2nd  readinq  o-f  Channel  0  sent  to  comouter 
;  Cn'J  REN  «■«  This  second  readinq  oF  Channel  0  is  indicated  bv  '  CAO ' 

'C~'’  REN  Third  readinq  would  be  CB::.etc.  imry*ir> 

:C30  LINE=LlNErl 

.C'^0  PRINT#  2  :  C0.C1.CA0,C2.C2.CA2.TIME 

;3<'0  FO-NOMF  (0)  /  10.  24#  •:C0+li;i'7)  I?  FAOwNQIhFCO) /10.24#  (CA0+10''7)  '  Converts  counts 

r  readings  to  >!tal  -freqs 

I3i;i  FisNOMF  ( 1 ) /'10-24#  (Cl-fl0-'7)  i*  F3=NQMF  (3) /iO.  24*  (C3+10''7)  '  Converts  to  trsa 
1320  F2-MQMF '2) / 10.24* (C2+10''7)  ®  FA2-NQMF (2) / 10. 24# ( C2* 1 0  7 )  '  Convers  to  frac 
1330  FIOLD-Fl  I?  F2QLDaF2  '?  F30LD=F3  FOOLO*FO 

1340  FAOOLO^FA  ( 0 i?  FA2CLDwFA<2)  Current  F  value  becomes  'old'  value 
1350  TEMPCQNST=93  I  Valid  only  -for  Colorado  Crystal  a/n  1  ###««« 

1360  REM  <*#  FOLLOWING  CONVERSIONS  NECESSARY  TO  MAKE  USE  OF  A  PREVIOUS  PRCGRAM 
1370  BMODEl^FO  !  1st  B— MODE  readinq  -for  crystal  under  test 

1330  TMaDE«F2  '?  FMQDE-F3  '  Used  to  give  TAPPROX  from  BS<C  MOCE  Obaratior  i?CHAM2-  3 
PRINT  CHR3  (27).r<'‘l  IL" 

1 39'.'i  3MQDE2*FA0  I  2nd  B'-MQDE  readinq  crystal  under  test 

1400  CMCDE»F1  I  In  between  C-MODE  readinq  for  crystal  under  test 

1410  TAPPRQXa-< (TM0DE-FMQDE-4700OO) /140)+TEMPC0NST  !  For  use  with  Colorado  s/n: 

1430  'CMQDEl'^BMQOEO /'2“CM0DE—470000  !  B— C  difference  freq.with  offset 

i  iT'!'  W=CflQDE-S000000  !  Deviation  in  Hz  tar  C-Mode 

144‘.)  P=R 1 .  D-*-2)  I  Inserts  A<0)  value 

1450  FOR  J=1  rC  D 

i4i3.i  F"“P-‘R  ( J  +  l .  D-^2)  tX"  J  '  Forms  polynomial  appro;:  i  mati  on  value 
14’;i  NETT  J 
1  430  CDEV=CMCDE-COLD 
1490  COLDaCMOOE 

:  5'.''  '  ;f  aBS  (CDEV)  30  THEN  BEEP  10.200 

1510  IF  ABS  (BM0DE1-BM0DE2)  300  THEN  BEEP  99,200  ELSE  1320 

1520  PRINT  USING  1330  :  L INE.  BMCDEl ,  BMODEl -BM0DE2.  CMOCE.  CDEV .  TuP=  POi  .  ’  I  ME  .  ' 

'5.  lo  o 

133'.;  image  4D.  1  X  .  3D.  3D.  IX.  S2D.  3D.  1 X ,  SD.  3D.  1 X  .  32D.  3D.  1  X  .  S3D.  2D.5X.5C.  4<.7^r. 

1 547:  PEM  #######««. 

ISaO  REM  *  THE  FOLLOWING  SUBROUTINE  WILL  PLOT  THE  INDIVIDUAL  =QrjTS  uT'  =22’’  « 

Sfl  REM  XXXXtXttXXXXXXXXXXXXXSXXXXXXXXXXXXXXXtXtXtXxxtvxxtXTtXttYtrXXXXXXvtxx 

I5”'i  (W-P- '  *20  '  xxxval  id  only  for  3KALE*3  U  “reqaj  nh; 

1530  MOVE  X.V 

.5®'i  IF  ■OLD-;.  X  THEN  1600  ELSE  l6lO 
I  j'.'O  LABEL  SEEP  10,100  i?  GOTO  1620 

111"  LABEL  "o"  '?  SEEP  13,100 
li20  XCLDwX 

:=3:’  R£M  «x**  A  long  WAIT  can  Pe  inserted  here  to  slaw  acwn  orintout  --ata 
1j4'.!  cm  lEVi*  1  GOTO  2080  '  Used  to  DUMP  GRAPHICS 
l.-^5'.'  GOTO  1630  '  Idle  loop  for  cvcle  timer 

1  s  af;  REN  KxtXXXxxsxXXXXXXX.tXXXXXXXXXXXXXXXXXtXXXXXtXlXXttttoitXtxtTXx'rxxt'c;** 

luT.-i  REM  t  THE  FOLLOWING  WILL  PLOT  THE  FRAME  AND  THE  AXES  i 

163"  REN  < XX** •t**X********X**X********X*X***** ********* **************** 

1;=0  CLEAR  ■}  GCLEAR 
■  -O':  =’L0TTSR  IS  1 
I'l'.'  il  r^LE-s^ 


C-3 


fr. 


‘4.:^ 


:  ~'30 
1  ~-o 


1370 
1 34';' 
l  530 
1  960 
1370 

isao 

:  3’’o 

'.910 
1920 
1  9T0 
1940 
I  930 
1 9t'.') 
1  ■  '  .1 
:  4  3i;i 
:  9<?o 


-■  .  '.I 
.'4<;i 


a>.i 


oaiid  +ar  Colorado  CrvstJl  S/'o  l;madi- 


5ri-;i--HALL 

L-JCATE  j'l' .  2-j':  .  3’.- ,  a~ 
i^R-rtriE 
utG 

<:'iN/'= r900';i  '?  xmin=-3  ;i(;jo 

CSIZS  4 

3CALE  XilAX.  aMIM,  -3KALE.  3KALE 
fAXia  XMIN. . 3. -SKALE, 3KALE 
LAXES  COO. 1 . XMAX.-3KAUE.3. .3.4 
nC'v'E  xMAX-aOOO.  -  (  l .  43*SKALE) 

csi:e  4 

LABEL  ‘OFFSET  T-F  FREQUENCY ':HZ ) 

MOVE  XMAX- 1 4000. - ( . 3*SKALE) 

'C3IZE  4 

LABEL  'RAW  DATA;  " ; DTE* 

MQ'OE  XMAX-aoOO.  1.37*SkALE 
CSI3E  4 

LABEL  "APFRax  TEMPERATURE  <C) " 

MQ'YE  XMAX  +  1300.-(  .  93*SKAL£) 

OES 

LDIR  90 

LABEL  "-'F/F  FREQ  DEVIATION  (10 
MOVE  XMAX+1600. .93»SKALE 
LABEL  "-a" 

MO'.'E  XMAX.O 
LINE  TYPE  3 
DRAW  XMIN.O 
r»-33 

FOR  1=0  TO  31  STEP  1 

X  =  X^1AX■^30-I  ^  XMAX-XMIN)  /31 

MOVE  X.1.03*SKALE 

label  r*31il 

MOVE  X..93*3KALE 

LABEL 

NEXT  I 

RETURN  '  ftet'jrn  to  1030  attar  plotting  CRT  frama 
printer  is  703  a  PRINT  ij  PRINT  i?  DUMP  GRAPHICS 

GOTO  l(.i3o  '  Loop  back  to  normal  data  logging  after  DUMP  GFAPhI''- 
3NB 


) 


C-4 


5  MISSION 

\  of 

^  Rome  Air  Development  Center 

n  HAVC  plum  and  vczaUti  A.&6&aA.ch,  d&v&lopm&nt,  tzM  and 
»  ^ej-zclad  cx.c.qiuj,itioYi  pKaq/uum  In  i^apponX  ot$  Command,  Contxoi 
2  Cotmuju.cadN.om  and  InteltCgcnce.  (C^T)  actcvlZcz^.  Technical 
%  and  cnqlmcaing  6uppo/U  mithin  oAzoi  o/^  technical  competence 
S  ij>  provided  to  ESV  Pfiogaam  O^^-ices  (PeJ-i)  and  othen.  ESD 
2  element^.  The  principal  technical  mli4lon  ateoi  a/ie 
L  comrnanicatcom ,  electkomagnetic  guidance  and  control,  ia-t- 
f  velUance  o(i  gaound  and  aenoipace  objects,  intelligence  data 
»  collection  and  handling,  information  6y4,tein  technology, 

2  ionospheric  propagation,  solid  state  sciences ,  microwave 
^  physics  and  electronic  reliability,  maintainabilitu  and 
s  compatibility. 


END 

FILMED 

3-85 

DTIC 


